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Abstract. We investigate certain measures induced by families of non- 
intersecting paths in domino tilings of the Aztec diamond, rhombus tilings 
of an abc-hexagon, a dimer model on a cylindrical brick lattice and a growth 
model. The measures obtained, e.g. the Krawtchouk and Hahn ensembles, 
have the same structure as the eigenvalue measures in random matrix theory 
like GUE, which can in fact can be obtained from non-intersecting Brownian 
motions. The derivations of the measures are based on the Karlin-McGregor 
or Lindstrom-Gessel-Viennot method. We use the measures to show some 
asymptotic results for the models. 



1. Introduction 

We begin by summarizing some facts about the Gaussian Unitary Ensemble of 
random hermitian matrices. We will see analogues of these in the random tiling 
problems discussed below. The Gaussian Unitary Ensemble (GUE) is the proba- 
bility measure 

(1.1) J_e-t^A^'dM 

on the space of all x hermitian matrices, which is isomorphic to M^^, and 
dM is the Lebesgue measure on this space. The induced measure on the N real 
eigenvalues Ai, . . . , Aa? of M is given by, [S^ , 

1 ^ 

(1.2) </>^,GUE (A)rf^A = — Ajv(A)^ TT e-^Wx 
where 

(1.3) AA,(A)=det(Af^'=),Vi= n 

l<j<k<N 

The probability of finding m eigenvalues in infinitesimal intervals dXi,...,dXm 
around Ai, . . . , A„i is given by Rm.NiXi, . . . , Xm)dXi, . . . , dXm, where Rm,N is the 
m-point correlation function 

Nl f 

(1.4) R,n,N{Xl, ■ ■ ■ , X,n) — jrri rr / 07v,GUE (Ai , . . . , Ajv)dA„j+i , . . . , dAw. 

[N -my. Jriv-„. 
The correlation functions are given by determinants, 

(1.5) Rm,N{Xu ■■■,Xra)= det(i^w(A„ Aj))™=i, 

where 

/-, ( \ KN^ihN{x)hN-l{y)~hN~l{x)hN{y) f 2\i/2 

(1.6) KN(x,y)^ e 

kn X — y V / 
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and hn(x) = K„a;" + . . . are the orthonornial Hermite polynomials. 

We can think of the eigenvalues as points on R. Using the formulas above and 
asymptotics for Hermite polynomials we can obtain a limiting random point process 
with correlation functions 

^sin7r(a;i 



(1.7) 



i?,„(xi, . . .,Xra) = det(- 



r(x2 "^j) 



3 I \rn 
)i,j=l- 



Hence we obtain a determinantal random point process given by the sine kernel, 
|62[ . To get this limit we rescale so that the mean distance between the eigenvalues 
(points) equals 1 and then use, 



sin 7r(^ — 77) 




V2Np{u) 



7r(^ - v) 



where p{u) = ^-^7(4 — u^)^ is the Wigner semicircle law which describes the 
asymptotic density of the eigenvalues. 

Let us recall two limit theorems for the fluctuations of the eigenvalues. Denote 
by #(u,w) the number of eigenvalues in the interval [u^N/2, v^jN/2], u < v, 
\u\,\v\ < 2. Then, 



(1.9) 



converges in distribution to the standard normal, []62[ . 
— maxi<j<7v Xj is the largest eigenvalue, then 



If A, 



(1.10) I 

as — * 00, ^ £ R, where 
(1.11) 



■AT 



A. 



- V2N 



F{0 



F{0=det{I^A)L2^^^oc)- 



Here A is the operator on L^(^, 00) with kernel (the Airy kernel) 



(1.12) 



A{x,y) ^ 



Ai (a:)Ai'(y)-Ai'(x)Ai jy) 
x-y 



p6[. The distribution function (1.11) is called the Tracy- Widom distribution. It 



follows from (1.5) and the Fredholm expansion that 



(1.13) 



E 



N 



N 



11(1 + 



N ^ 

V— / det{KN{x,,Xj)g{xj))lj^^d''x 



det(/ + KNg)Ly 



If we take g{x) — ~X{t,oo)ix) in (1.13) and use the asymptotics of the Hermite 
polynomials close to the largest zero, we can prove (1.10). 



The GUE eigenvalue measure, (1.2), was obtained above from the measure (1.1) 
on random hermitian matrices. We will now show how we can obtain (1.2) in a 
completely different way using non-intersecting Brownian motions. This type of 
problem has been studied under the name of vicious walkers or domain walls in the 
statistical physics litterature, see ||l^, |19 , |Q and references in these papers. 
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Consider N 1-dimensional Brownian motions starting at the points 0, . . . , — 1 at 
time and ending at the same points 0, . . . , N — 1 at time 2T. Let pn.t{xi, . . . , x^) 
denote the probabiUty density that at time T the particles are at the positions 
xi < ■ ■ ■ < xjv conditioned on the event that the paths have not intersected in the 
whole time interval [0, 2T]. If pt(a;, y) = (27rt)~-'^/^ exp(— (x— ?/)^/2t) is the transition 
kernel for Brownian motion, then, by a theorem of Karlin and McGregor, [Q, see 
also @, 

(1.14) pm,t{xi,. . . ,xn) = -^det{pT{j - l,Xk))f,k=i^(^^{PT{xj,k - 

ZjN 



^{det{pTij-l,Xk))f^,__ 



where 



(1.15) ^^ = ^f {det{pT{j-l,Xk))f,,^,)'d''x. 

Note that, because of symmetry, we can consider ( |1.14[ ) as a probability measure on 
and remove the A^! in (1.15). It follows from |67|| that the probability measure 
on with density (1.14) has determinantal correlation functions analogous to 
( p^ ) but with a different kernel, se e pof . We can now obtain GUE as follows. Let 
T — > CX3 and rescale the Xj'.s by vT so that they do not move away to infinity. 
Then, 

(1.16) lim PN.rixiVf, xnVt) (j)N.GVE (x). 

T — >-oc 

To show this note that, 
(1.17) 



det(pH, - l,x,))f,,^, = ^^-^ n e«^-i)^+^?)/^^det(e(-i)-/2^),^, 



and use the formula for a Vandermonde determinant. The choice of initial and final 
positions made above is not necessary for the result but simplifies the computations. 
Thus, we need not look upon (^]^) as something which necessarily comes from 
random matrices. For another relation between random matrices and Brownian 
motion see p^ . 

The present paper can be seen as a continuation of the paper where several 
examples of analogues of GUE on a discrete space was given. In these ensembles 
we have analogues of the results for GUE discussed above but we obtain the so 
called discrete sine kernel, in the limit instead of the ordinary sine kernel. The 
common theme of the present paper is non-intersecting paths, which are discrete 
analogues of the non-intersecting Brownian motions just mentioned, and thus it 
is reasonable to expect that we will have features analogous to those of GUE. In 
section 2 we will consider random tilings of the Aztec diamond of size n, psf , which 
can be described using certain non-intersecting paths. By using so called zig-zag 
paths in the tiling, we obtain the Krawtchouk ensemble, which can be used 
to analyze several properties of the random tiling. By a result in Jsof the shape of 
the so called temperate region in a random tiling is closely related to the corner 
growth model in , which is a generalization of the longest increasing subsequence 
problem for random permutations. This gives a new approach to the asymptotic 
fluctuation results in B and involving the Tracy- Widom distribution. We will 
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also study the fluctuations of the domino height function which describes the dihng. 
It has Gaussian fluctuations with variance of order log N, a fact that is related to 
the Gaussian fluctuations in the number of eigenvalues in an interval in a GUE 
matrix. This type of result has been conjectured in |5^. We will also discuss 
the relation between the equilibrium measure for the Krawtchouk ensemble and 
the arctic ellipse. The corner growth model can be generalized further and this 
leads to the so called Schur measure introduced on Q. The Schur measure can 
also be analyzed using non-intersecting paths as will be demonstrated in section 3. 
We will also, in section 4, consider rhombus tilings of a hexagon, which are 
related to boxed plane partitions. These tilings can also be described by certain 
non-intersecting random walk paths and the intersection of these paths with a fixed 
line leads to the so called Hahn ensemble. In this problem we will not compute 
the detailed asymptotics, but we will discuss the equilibrium measure for the Hahn 
ensemble and its relation to the arctic ellipse phenomenon. Finally, in section 5, 
we will analyze certain aspects of a dimer model on a brick (hexagonal) lattice on a 
cylinder. Here we also have non-intersecting paths but the number of paths is not 
fixed like in the other examples. The methods used are very close to the arguments 
used to compute the correlation functions in [6^ 

There are many papers in the statistical physics litterature related to the present 
paper, e.g. §,]l§,@,|l|,|§,|§,@,g|,|6§and0 Connections between 
random permutations and the so called random turns model, which gives certain 
non-intersecting paths has beeen discussed in [Q, and |^. Other relevant 
papers are [||, ||, |8|, || and [g. 



2. The Aztec diamond 

2.1. Basic definitions. The Aztec diamond, An, of order n is the union of all 
lattice squares [m, to-|- 1] x [l,l + l],m,l £ Z, that he inside the region {{x,y) ; \x\ + 
\y\ n + l}. A domino is a closed 1 x 2 or 2 x 1 rectangle in with corners in Z^, 
and a tiling of a region i? C by dominoes is a set of dominoes whose interiors 
are disjoint and whose union is R. Let T(A„) denote the set of all domino tilings 
of the Aztec diamond. 

We can equivalently think of a tiling as a dimer configuration. Consider the 
graph G with vertices at {^,^) + Z'^ and edges between nearest neighbour vertices. 
A dimer is simply an edge in G, and if the edge goes between the verices vi and 
V2 we say that the dimer covers vi and V2- Let Gn be the subgraph of G where all 
vertices lie in An- A dimer configuration in Gn is a set of dimers in Gn such that 
all vertices are covered by exactly one dimer. This is clearly equivalent to a tiling 
of An via the identification: a dimer between vi and V2 corresponds to a domino 
covering the two lattice squares with centers vi and V2- 

Colour the Aztec diamond in a checkerboard fashion so that the leftmost square 
in each row in the top half is white. A horizontal domino is north-going (N) if 
its leftmost square is white, otherwise it is south-going (S). Similarly, a vertical 
domino is west-going (W) if its upper square is white, otherwise it is east-going 
(E). Two dominoes are adjacent if they share an edge, and a domino is adjacent 
to the boundary if it shares an edge with the bundary of the Aztec diamond. The 
north polar region is defined to be the union of those north-going dominoes that are 
connected to the boundary by a sequence of adjacent north-going dominoes. The 
south, west and east polar regions are defined analogously. In this way a domino 
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tiling partitions the Aztec diamond into four polar regions, where we have a regular 
brick wall pattern, and a fifth central region, the temperate zone, where the tiling 
pattern is irregular. 

We will now define a one-to-one mapping from T(j4„) to families of n non- 
intersecting lattice paths. Consider an S-domino which we place with corners at 
(0,0), (2,0), (2, 1), (0, 1). Draw a straight line from (0,1/2) to (2,1/2). In this way 
we get a piece of a path, and we do this for all S-dominoes. Similarly we can put 
a W-domino so that it has corners at (0, 0), (1, 0), (1, 2), (0, 2), and then draw a 
straight line segment from (0,1/2) to (1,3/2). Finally, on an E-domino placed at 
the same position we draw a straight line from (0,3/2) to (1, 1.2). Do this for all 
W- and E-dominoes. We do not draw any line on an N-domino. Given a domino 
tiling of An we draw lines on the dominoes as just described. We claim that this 
gives n non-intersecting paths starting at Aj = (n -I- 1 — j, i — j) and ending at 
Ej = (— n — 1 < j <n. We call these paths DR-paths of type I (after D. 

Randall, |Q, p. 277). To prove the claim we argue as follows. Consider the black 
lattice square to the right of Ej. It has to be covered by a W- or am S-domino. In 
both cases a path will end at Ej From the checkerboard colouring we see that we 
nust obtain connected paths. Similarly, if we consider the white lattice square to 
the left of Aj it can only be covered by an E- or an S-domino. Hence a path must 
start at the point Aj. Clearly, by construction, the paths are non-intersecting and 
the claim is proved. 

A convenient coordinate system for describing the paths is what we call co- 
ordinate system I (CS-I). As origin we take {n + 1,1/2) and as basis vectors 
6/ = (—1,-1), f/ = (—1,1). Let £/ be the integer lattice in CS-I. The type I 
DR-paths are walks in £/. They take steps (1,0), (0,1) or (1,1) and they have 
starting points (fc,0) and endpoints (ji + 1, k), 1 < k < n. Thus we obtain a map 
from domino tilings of An to families of n nonintersecting type I DR-paths in £/ 
with the specified initial and final positions. This map is a bijection. To see this, 
fill in with dominoes along the paths using the marked tiles. This is possible since 
the paths do not intersect. If we have a white lattice square that is covered by a 
domino, then the black square to the right is also empty, since otherwise the paths 
would not be connected. Similarly, if a black square is not covered by a domino, 
the white square to the left is not covered either. Hence, the squares that are not 
already covered can be covered by N-dominoes. Clearly this gives an inverse. 

We can also define type II DR-paths which are complementary to the type I 
paths. In this case the S-dominoes are unmarked, whereas the N-dominoes have a 
horizontal segment in the middle. Furthermore we interchange the marking on the 
W- and E-dominoes. In this way we obtain paths from Aj = (— n — 1 + j, i — 1/2) 
to Ej — {n -\- 1 — j,j — 1/2), 1 < j < n. In coordinate system II (CS-II) which has 
origin (— n — 1,-1/2) and basis vectors e// = (1,1), f// = (1,-1), we obtain the 
same type of lattice paths as before. 

Call the top type I DR-path the level- 1 path. It is clear from the definitions 
above that the north polar zone is exactly the part of the Aztec diamond above the 
level-1 path, i.e. all dominoes, which have to be N-domonioes, that lie above this 
path. 

Let T e T{An) be a tiling of the Aztec diamond and let v{t) denote the number 
of vertical dominoes in t. We define a probability measure on T{An) by letting the 
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horizontal dominoes have weight 1 and the vertical dominoes weight w. Thus, 

(2.1) PM = ^. 

If we take ?i; = 1 we obtain the uniform distribution on T(A„). This can alterna- 
tively be viewed as a probability measure on the DR-paths, where we put the weight 
1 on the steps (1, 1) (in CS-I or CS-II), which correspond to horizontal dominoes, 
and the weight w on the steps (1,0) or (0,1), which correspond to vertical dominoes. 
The weight of n given non-intersecting DR-paths is the product of the weights on 
all steps and equals w'"^'^\ if r is the tiling defined by the paths. The weight of a 
set of non-intersecting DR-paths is the sum of the weights of all the elements in 
the set. 

Next, we will define another type of paths, the so called zig-zag paths, [|l5j, in 
the Aztec diamond. Consider the sequence of white squares with opposite corners 
= (— r + k,n + I ^ k — r), fc = 0, ...,n-|-l, where r, 1 < r < n, is fixed. A 
zig-zag path Zr in An is a path of edges going around these white squares. When 
going from Q^. to Ql._^_i we can go either first one step east and then one step south 
(an ES-step), or first one step south and then one step east (an SE-step). A domino 
tiling T e T{An) defines a unique zig-zag path Zrij) from Qq to Qn+i if we require 
that the zig-zag path does not intersect the dominoes. There will be exactly r 
ES-steps, and hence n + 1 — r SE-steps along the zig-zag path. This can be proved 
using the domino height function defined below. 

We associate the point (r, n — fc) in CS-I with the step Suppose that we 

have ES-steps at the points {r,hj), 1 < j < r, in CS-I, i.e. Q^_i^.Qn_fi._^_i are ES- 
steps. Then the zig-zag path is mapped one-to-one to {hi, . . . , hr), where < hi < 
■ ■ ■ < hr < n. This specifies the zig-zag (particle) configuration (hi, . . . , hr); we 
write p{Zr) — {hi, . . . , hr)- We can also associate the step Q^Qk+i ^i^h the point 
{n-\-l — r,k) in CS-II, and we will then have SE-steps at the points {n-\-l — r,n — kj), 
^l^jl^n+l — r, where ki < ■ ■ ■ < kn+i-r- We call (fci, . . . , fc„+i_r) the zig zag 
(hole) configuration, and write h{Zr) — {ki, . . . , kn+i~r)- The next lemma gives 
the relation between the DR-paths and the zig-zag paths. 

Lemma 2.1. The points {r, hj) in CS-I are the last positions on xj — r of the type 
I DR-paths starting at {k, 0), 1 < k < r in CS-I. Similarly, {n + I — r,n — kj) are 
the last positions on xji ^ n + I — r in CS-II of the type II DR-paths starting at 
{k,0), l<k <n-\-l-r. Also, 

(2.2) {hi, ...,hr}U {ki, kn+l^r} = {0, . . . , u} . 

Proof. If we have an ES-step around a white square, then this square is covered 
by an S- or a W-domino. In both cases it follows, from the definition of the type 
I DR-paths, which are walks in the integer lattice Cj in CS-I, that the DR-path 
must intersect the S-step in the ES-step, and after that go to a point with a larger 
^/-coordinate. Similarly, if we have an SE-step around a white square, then this 
square is covered by an N- or an E-domino, in which case a type I DR-path does 
not intersect neither the S- nor the E-step. The proof of the second statement in 
the lemma is analogous, and (|]^) follows from the definition of the particle and 
hole configurations, and the definition of CS-I and CS-II. □ 

If p[Zr(T)] = {hi, . . . ,h^.), 1 < r < n, then the position of the rightmost particle, 
h!^, describes the level-1 type I DR-path. As noted above, the region above this 
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DR-path is the north polar zone, and hence we can investigate the shape of the 
north polar zone using the positions of the rightmost particles. We will return to 
this in sects. 2.3 and 2.4 below. 

Let us recall the definition of the (domino) height function associated with a 
given tiling, Let u and v be two adjacent lattice points (vertices), in the basic 
coordinate system, such that the edge connecting them is not covered by a domino. 
If the edge from u to v has a black square to its left, h{v) — h{u) + 1, and if it 
has a white square to its left, h{v) = h(u) — 1. Note that the value of the height 
function is uniquely determined apart from an overall additive constant. We can 
fix it by requiring h(n, 0) = 0. If w and v are two adjacent lattice points, then 
\h{u) — h{v)\ = 3 if the edge is covered by a domino, otherwise \h{u) — h{v)\ = 1. 
From this it follows that h{Ql) = 2n - (2r - 1), /i(Q;+i) = 2r - 1 and 



(2.3) MQD-MQfe+i) 



-2, if ES-step 
2, if SE-step. 



Consequently, we can use the zig-zag configurations to determine the height at a 
given point. From ( ^.3| ) it follows that there are exactly r ES-steps in Zrir). 

2.2. The Krawtckouk ensemble. The Krawtchouk ensemble, |Q, is a probabil- 
ity measure on {0, ... , K}^ defined by 

1 ^ /K\ 

(2.4) FKr,N,KAh] = ^ A^(/i) n , 

where 0<p<l, q = l-p, l<N<K,h= (/ii, . . . , /iw) € {0,...,K}'^ and 

(2.5) Z^,K,p = iV! ( n I Kl^ipqr^''-'^/'. 

Set w{x) — < X < K and let {pkix)}^^Q be the normalized orthogo- 

nal polynomials with respect to w{x) on {0, . . . , K}, i.e. 

K 

(2.6) y^^pj{x)pk{x)w{x) = 6jk- 

These are multiples of the ordinary Krawtchouk polynomials, [ p4[ , and have the 
integral representation 

^^"^''^„„^-«/2 1 f (l + qzYil-pz)^-'^ dz 



(2.7) ^ J (M)-"-^ 

where 7 is a circle centered at the origin with radius < min(l/p, 1/q). The measure 
(^^) has determinantal correlation functions, |52|, |]67|| , 

(2.8) det{KKr,N,K,p{^^' 

where the Krawtchouk kernel is given by 

(2.9) KKr,N,K,p{^,y) = ^Pn{x)pn{y)(w{x)w{y)y^'^ 

n=0 



KN^lPN{x)pN-l{y) - PN-l{x)pN{y) , I . , 
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1/2 

where k„ — (p<z) "^^ is the leading coefficient in Pnix). 

We win now prove, using the DR-paths. that the measure on the zig-zag config- 
urations induced by the probabihty measure (2.1) on the tilings is the Krawtchouk 
ensemble. In the special case of uniform distribution on the set of tilings this was 
proved in using formulas from [T5|] . 

Theorem 2.2. Fix r, 1 < r < n, and let h = {hi, . . . , hj,), where < hi < ■ ■ ■ < 
hr < n be given. Then, 

(2.10) P[p{Zr{T)) ^h]^ rlPKr,N,KAh], 

where q — 'w'^{l + w'^)~^ . Hence, if we disregard the ordering of the particles in the 
zig-zag particle configuration, the probability of h is exactly Ppcr.N.K.qih]. 



Proof. By lemma 2.1 we have type I DR-paths from (r + 1 — j, 0) to {r,hj), 
1 < i < in CS-I, and type II DR-paths from (j, 0) to {n + 1 ~ r,n ~ kj), 
l<j<n + l — rin CS-II, where fci < • • • < kn+i-r and ( pT^ ) holds. Together 
these describe the whole domino tiling. Let u![h] be the weight of all the type I DR- 
paths between the specified points and u' [h] the weight of all the type II DR-paths 
between the given points. Then, 

(2.11, mzArfl-Hi^ "'"'"'"■ii,, 

The quantities LL}[h] and Lu'[h] can be computed using the Lindstrom-Gessel-Viennot 
method, [|ll, Hi, see also which is a development of the Karlin-McGregor 
result in a discrete setting. We want to compute the weight of a path that takes n 
steps to the right and m steps up. Let a be the number of (1, 0) steps, b the number 
of (0, 1) steps and c the number of (1, 1) steps. Then, n ~ a + c and m = b + c. 
The number of paths with a given number of steps a,b,c equals 

{a + b + c)\ {n + m — c)l 

alblcl {n — cy.{m — c)\c\ 

Note that c can take all values between and min(n, m). The total weight of all 
possible paths from (0, 0) to {n, m) is thus 

min(n,m) 

(2.12) w{n,m)= ^ -y^--^l—2_y^ 



^ (n + TO-c)! ^„„+„,„2e 

(n — c)!(m — c)!c! 

c— ^ ' ^ ' 



If we use the Pochhammer symbol (a)k — ci{a -I- 1) . . . (a + fc — 1), (a)o = 1, then 
(2.13) w{n,m) — j — > [m — c + l)m{n — c+ 1) 



n\ ^ — ' c! 

The Lindstrom-Gessel-Viennot method now shows that the total wight of all pos- 
sible r non- intersecting paths from (r -f 1 — j, 0) to (r, hj),\ < f <.r, is 

(2.14) ^\h\=A^X{w{i-\,h^))l^^i. 

Similarly, the total weight of all possible n + l — r non-intersecting paths from {j, 0) 

to (n + 1 — r, n — kj) is 

uj'[h] = det(w(n + l-r-i,n- 

= det(w{i-l,kn+2-r-j))ljlr- 
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If we set r' — n + 1 — h'j = n — kn+2-r-j, then 



(2.15) 



c^'H =dct(w(z-l,/i;.))[,,=i, 



which has exactly the same form as ( 2.14 ). This is the advantage of using both 
types of DR-paths. 
Now, 



(2.16) 



det 



To see this, insert ( ^.13 ) into the left hand side of ( ^.16| ), and use the muhihnearity 
of the determinant to obtain 

,,j-l+hj 



IIt— TV E n(*-^Oc.^det((/.,-c. + i)._or,., 



which equals the right hand side of (2.16) since. 



c=0 c=0 ^ ' 

From ( p.l5D we obtain, after some manipulation, 



n+l — 



'■I'l = I in iJi^ n 



W 



n—kj 



i=i 



Lemma 2.2 in |38| shows that if si < • • • < sjv and ri < ■ ■ ■ < vm and the union of 
these two sets of numbers is exactly {0, . . . ,N + M — 1}, then 



' N+M-l 



N 



(2.17) AM(r) ^ n ^-11 



1 



i=i 

If we use this formula, we obtain 



l^s,liN + M-l-s,)\ 



uj[h]Lj'[h] = u;-'^('-i)(l + u;2^»("+i)/2-"'-+'-('- 



,,2/i,- 



By this formula and (2_J) we obtain 



U^ihWih] = (l+^„2)„(„+l)/2 - ^^(/^)2-Q > 

where g = w^(l + w^)^^. Hence 

(2.18) t^[;i]cj'[;i] = (1 +ii;2)"(»+i)/2, 

0</ii<-</i,.<n 



and the theorem follows from (2.11). 



□ 
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As a corollary we obtain the following result first proved in ||T5| , see also [|6) . 

Corollary 2.3. The number of elements in T(A„) is 2"("+i)/2^ and the probability 
of having 2k vertical tiles, < k < n{n + l)/2, is 

The probability of having an odd number of vertical tiles is zero. 



Proof. The first result follows by putting w — 1 in (2.18), and (2.19) follows by 



expanding the right hand side of (2.18) using the binomial theorem. □ 



2.3. Asymptotics in the Krawtchouk ensemble. There is an equilibrium mea- 
sure associated with the Krawtchouk ensemble, see sect. 2.2. in |3^, and sect. 
4.2 below for some more details. If r — > cx) and n — > cx) in such a way that 
r/n —> t G (0, 1), then the expectation of the discrete measure ^ Sj^=i ^hj/n con- 
verges weakly to the equilibrium measure ut^q{x)dx, i.e. the equilibrium measure 
gives the asymptotic distribution of the particles. The equilibrium distribution is 
scaled so that its support is contained in [0, 1]. The equilibrium measure also gives 
the asymptotic distribution of the zeroes of the Krawtchouk polynomials scaled to 
[0, 1]. See for this result and explicit formulas for the equilibrium measure and 
its support. Since the position of the rightmost particle determines the boundary 
of the north polar zone (and we can make an analogous analysis for the other polar 
zones or use symmetry), we can use the equilibrium measure to prove the arctic 
ellipse theorem, and jl^. The arctic ellipse theorem of Jockush, Propp and 
Shor says that the boundary of the temperate zone, scaled by l/n, converges in 
probability to an ellipse as n — > oo. 

Theorem 2.4. // we scale the Aztec diamond by l/n in the original coordinate 
system, then the boundary dTn of the temperate zone in a rescaled andom tiling 
of An under the probability measure ^2.^ ), converges in probability as n ^ oo , 
r/n t Cz (0, 1), to the ellipse E, 

2 9 
X V 

(2.20) — + — = 1, 

p q 

in the sense that P[dist{dTn, E) > e] for any fixed e > 0. Let dist i{dTn, E) be 
the maximal distance from a point on dTn inside E to E, and dist o{dTn, E) be the 
same thing but from a point outside E. Given e > 0, there are positive constants 
/(e) and J(e) such that 

VaasnY>\\ogT[disti{dTn,E) > e] < -/(e) 

n — *oo n 

and 

limsup-logP[dtsto{dTn,E) > e] < -J(e) 

n — *oo n 

Proof. We will indicate how the shape ( ^.2Cl| ) is obtained. The large deviation 
formulas, which imply the convergence in probability to the ellipse, follow from 



theorem 2.2 in |37 . See sect. 4.2 below for some more details in the analogous 
result for rhombus tilings of a hexagon. Let {xi,yi) be coordinates in CS-I and 
(a;, y) coordinates in the original coordinate system. Then, x = n + l/2 — xi — yi, 



y — l/2~ XI + yi. It follows from 1 14 that the right endpoint of the support of the 
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equilibrium measure is given by (3{t, q) — tp+{l — t)q + 2y/t{l — t)pq if < t < 1 — g. 
The connection between the position of the rightmost particle and the boundary 
of the north polar zone described above imply that the limiting boundary of the 
north polar zone must be the curve (0, 1 — g] 3 t {t, P{t, q)) = {xi,yi) in CS-I. 
Using the coordinate transformation we find that, in the original coordinate system, 
points on this curve satisfy (2.2C). □ 

Let i>[a,b]{h) be the number of particles in the interval [a, 6] in the particle 
configuration h. From the formula (2.3) and h{Q'^_^_^) = 2r — 1, we obtain 

n 

h{Ql) - (2r - 1) = Y.i'^iQ^j) - '^(Q.Vi)] 

j=k 

= -2u[0, n- k]+2{n- k+1 - u[0, n - k]) ^ 2(n - fc + 1) - 4//[0, n-k]. 
Consequently, 

(2.21) h{Ql) = 2{n-k + r) + l- 4:v[0, n - k]. 

Assume that k/n t and r/n ^t,0<T<ta,sn^ oo. Then, by (2.18) and the 
weak convergence of the particle distribution, 

h(O^) /-i— ^ 

(2.22) lim — = 2(l-r + i) -4i / utJx)dx. 

n^oo n Jq 

The precise form of the equilibrium measure is given in [ p^ and using this we can 
work out the asymptotic height function. We will not evaluate these integrals here. 
We can also obtain large deviation formulas (and estimates, compare lemma 4.1 
in |3^ ) for macroscopic deviations from the asymptotic (average) height function. 
This is analogous to the large deviation formulas for the Wigner semi-circle law, [Q . 
We will not develop the details, since it is very analogous to the corresponding ran- 
dom matrix results. See |10 for previous large deviation estimates and asymptotics 
for the height function. 

We will now analyze the fluctuations of the height function, or what is equiv- 
alent, by ( ^.2l| ), the fluctuations in the number of particles in an interval in the 
Krawtchouk ensemble ( ^!4| ). Let I = {b - L,b - L + I, ... ,b} C {0, K} he an 
"interval" of length L and let 1/(1) denote the number of particles in /, 

'^(I) = ; I <i < N and hi e I}, 

We want to prove the following result for the variance of i^(/), the number variance. 

Proposition 2.5. Assume that N/K t e (0,1/2], p = q = 1/2, b/K (3, 
[b — L) / K —t (3' < (3 as N^K,L —t oo, where (3 or j3' belongs to the interior, S, of 
the support of the equilibrium measure Ut.i/2- Then 

var{v{I)) 1 
log L TT^ 

where ^f3^p' is = I if both [3 and j3' belong to S and = 1/2 othe 



(2.23) 



This type of result was conjectured in ||5^. We will give the proof of this propo- 
sition below, which is rather long. The case p ^ q could also be worked out, but 
we stick with p = g = 1/2 for simplicity. Once we have this result we can apply the 
Costin/Lebowtiz/Soshnikov argument, ijll], |^l[], to prove that the fluctuations are 
normal. This type of results have been proved in other tiling models by Kenyon, 
see [g, @, 0. 
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Theorem 2.6. With the same assumptions as in the proposition, 

\/var[I) 

as K oo. i.e. we have convergence in distribution to a standard normal random 
variable. 

Proof. By a theorem in p. 8, see also the result ( ^.24 ) follows from (2.23) 
since we have a determinantal random point field, i.e. the correlation functions are 



given by determinants as in (2.8). The kernel ( p.9D defines a trace class operator 
(it has finite rank) /C, which satisfies < K. < I. This follows from (2.9) and 
the orthogonality (|2.6|). Hence the conditions in the Costin-Lcbowitz-Soshnikov 



theorem are satisfied. □ 



Combining this theorem with (2.21) we obtain the next theorem. 

Theorem 2.7. Take the uniform distribution on T(yl„) and < r < n/2 (the 
case n/2 < r < n is similar by symmetry) . Let = {—r + fc, n + 1 — fc — r) as 
before and let h(Ql,) be the value of the domino height function above this point. 
If r/n ^ t e (0,1), k/n n e (1/2 - ^i(l-t),l/2 + y^t^l - t)) = Ut or 
j /n — > k' G Ut, and \k — j\ — )■ oo, then 

(2.25) ^ ^ ^ ^(".1) 

4VC«,«'7r-2 1og|/c- j| 

Here ^k,k' = 1 if k,k' G Ut and ^k,k' — 1/2 if one of k or k' does not belong to Ut. 

Note that the condition on k, k' corresponds exactly to the condition that one of 
(3 or in proposition 2.5 belongs to the interior of the support of the equilibrium 
measure. 

We turn now to the proof of proposition 2.5 which is rather lengthy. The proof is 
based on sufficiently good asymptotic control of the Krawtchouk kernel, (2.9). For 
results on asymptotics of Krawtchouk polynomials see |3^. We state the needed 
result as a lemma, which we will prove later. Let p(^) be defined by ( p. 5^ ) below. 
We consider the case when p = 1/2 and N/K t. Then, p'(^) = Ut,i/2(C) is the 
equilibrium measure, which is suppported in [1/2 — \/t{\ — t), 1/2 + \/t{l — t)], see 



14 1 . A computation gives. 



TT ^l/4-t(l -t) 

Lemma 2.8. Consider the Krawtchouk kernel ^2.!^ with p = 1/2 and let 6 > 
(small). Sett = N/K and assume <t< 1/2. If\x/K~l/2\ < yJt{l~t)-5 and 
\y/K - 1/2| < yjt{\ - t) - 6, then th ere is a constant C , independent of N, K, x 
and y such that 

(2.27) 

1 , \^-y\ 



sm7rK{pix/K)-p{y/K)) 

{X - y)KKr,N,K,l/2(.X, y) 



< c 



Vk k 



Also, if \x/K — 1/2| < \Jt(\ — t) — 5, there is a constant C such that for y > x, 
(2.28) 

\ix ~ y)KKr,N,KM2{x.y)\ < CuAn{K''\ \t{l -t)- {y/K - l/2)\^l^)j^^^^. 
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There is also an analogous result for y < x. If 1/2 + \/t{l — t) < 1, we can obtain a 
much better estimate for y outside the support of the equilibrium measure; compare 
l ^.6di ) below. 

We will now prove proposition 2.5. 



Proof. Write IC{x,y) = -ft^Kr .jv,^, 1/2(2^1 2/)- It follows from ( |2.9| ) and the orthogo- 
nality (^.6[), that /C(x, y) is a reproducing kernel 



K 



(2.29) ^/C(*,j)/C(j,A:) =/C(*,A:). 



3=0 



Set / = {0, . . . , if} \ /. Then by the formulas (U|) and ( ^.29D , we see that 



jei \i=a / i.jei 



(2.30) =EE/C(*,jf = Ei + E2, 
where 

L K-b 

i=0 i=l 
h-L L 

Y.i = ^^]C{b- L-i,b- L + jf. 

i=0 j=l 

We will consider the case when f3 lies in the support S of the equilibrium measure. 
The contribution to Si comes from the right endpoint of the interval, and the 
contribution to E2 from the left endpoint. We will show that 

(2.31) Si = ^logi + 0(log(logL)). 

The s ame thing is true for S2, with an analogous proof, if /?' e S, so if we establish 

( 2.31 ) the proposition is proved. 

Let a{L) = [L/logL]. We spht Si into two parts 

L a(L) L K-b 

i=0 j=l 1=0 j=a{L) + l 

(2.32) = S'l + S'l'. 
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Now, by (|2.28D, 



L [(q-<5/2)A']-6 



1 



1 



1=0 j=a(L) + l 
oo [aK]-[Ki/^]-f, 

E E 

1=0 j = [(a-(5/2)A']-6 
CO [aK]-[K^'-^]+b 

E E 

i=0 j = [aA]-[A'i/2]_h 
oo K-b 

E E lifa- (6 + j)|l/2 (i+j)2' 

1=0 j=[aK]-[A:i/2]+b I V I V 



|ifa-(6 + j)|i/2 (j+j-)2 



where a = 1/2 + — t), the right endpoint of the support. The first sum is 

L oo 



< 



^E E 



1 



<Clog(logL). 



i=0 j=a(L) + l 

In the last three sums, we use YITLq + j)^ ^ l)i ^i^d it is then easy to 

see that the j-sums are < C. Thus 

(2.33) E'l' < C(l + log(logL)). 
We also split S'j^ into two sums 

a(L) a(L) L a(L) 

^'i = EE^(^-*'^+j')'+ E E^(^-*'^+j')' 

i=0 j=l i=a(L) + l j=l 

(2.34) =5i+^2. 



The second sum is estimated in the same way as S" using ( p.2q ), and this gives 
(2.35) S2 < C. 

To control Si we use ( p^ ), which gives 

sin2 nKipiib - i)/K) - p{{b + j)/K)) 



a(L)a(L) . 3 



i=0 j = l 
a{L) a(L) 



7r2(z+j)2 



(2.36) 



If we use sin^ a; = (1 + cos 2x)/2 and observe that 

a{L)a{L) 

(2-37) E E 7---yI = log^ + 0(log(logL)) 

j=0 j = l I.* 

we see that what remains to be proved is 

a(L) a(L) 

(2.38) E E 

1=0 j = l 



cos27rA'(p((6 - i)/i^) - p((6 + i)/ii:)) 



7r2(i + j)2 



< Clog(logL). 
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When we have this we note that (2.32) - (2.3S) imply (2.31) and we are done. 
To prove (2.38) we use summation by parts. Set S{0) — and 



S{j) = exp(-2^if zp(^)), l<j< a{L). 

r 

Consider the expression 
(2.39) 



a(L) 



a{L 



^ (m + 7')^ 

m=0 j = l ^ ■' ' 

whose real part is what we want. Summation by parts gives 



(2.40) 

a(L) 

V - 

^ (m + ]Y 
3=1 ^ ■" 



a(L)-l 



(m + a(L)): 



r5(«(i))+ E 



2(m + j) + 1 



(m + j)2(m + j + !)■ 



Since \S{a{ L)\ < - a{L), the first term in ( 2.4C ) gives a c ontri buti on < 1 to the 
expression ( |2.39| ). The contribution of the second term in ( 2.4C ) to ( 2.3£ ) is 



(2.41) 



< 



a{L) a(L)-l 

E E 



2{m + j) + 1 



=0 i^i + iY{m + i + 1)' 



:Sij). 



Fix an integer A > 1 and assume that A < j < a{L). Write j = kA + r, where 
< r < A, and fc < j/A. Then 



k A 



(2.42) 



Sij) = E E e-^'^^^^f^^) + E e-'^ 



Write ^ = {b + uA)/K. Then p( ''+"^+^ ) = p{C) + f p'(0 + 0{j^). Note that 
p' {C) > since we are inside the support of the equilibrium measure p' . Also, 
p'(0 < 1 by ( |2.26| ). Using this we see that 



(2.43) 
Now, 
(2.44) 



K 



A 



J2^-2.K 



ip(^)-2i:ivp' (^) 



since < p' (X) < 1. Similarly 

A 



(2.45) 



Ee-- 



<C, 



A3 



<C(1+-) 



Combining (|2l^) to (|2l5|) we obtain 



A3 1 A2 

|5(j)|<CMl + -^)<Cj(^ + — ) 
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for A < j < a{L). Since evidently 1 5*0)1 < j for 1 < j < A, we obtain the estimate 



a(L) a(L)-l 

E E 



2{m + j) + 1 



:S{j) 



<ClogA + C(^ + ^)logi 



of the expression in (2.41 ). If we take A — [logL], we see that the absolute value 
of the expression (|2.39|) is < log(logL), which gives the bound (|2.38|). □ 



We still have to prove lemma 2.8. 
Proof. Let p = q = 1/2 m (|2.9|) and make the change of variable z = 2w. This 



gives, 



Set 



Pn{z) = 



2ni 



(2.46) AN{x,y) 



KN 



-1/2 



K 

N ~ 1 



-1/2 



Gn{z;x) = {l + zY{l-z)^ ""w ^ and LNix,y) ^ {x~y)KKr^N,Ks/2{x,y). Then, 



LN{x,y) = AN{x,y) 



dz 



dw 



2'Kiz L, 2'Ki'w 



GN(z]x)GN{w;y){w - z). 



Set i = x/K and -q = y/K and define f{z) = Clog(l + z) + (1 - ^ log(l - z) - 
ilogz, which is the relevant function in the saddle point argument. The equation 
J'{z) = has the solutions z± = (1 - t)^i(^ - l/2±i^t{l-t) - (C - l/2f). Let 
Tc — \z^\ — \/t{l — t)~^ . We choose 7 to be the circle with radius rc, and write 



; exp(±z0c(O)j where 



(2.47) 



We assume that ^ > 1/2 (the other case is similar by symmetry). Also, we will first 
assume that ^ < 1/2 + i/t(l - t), i.e. we are inside the support of the equilibrium 
measure. It follows that < 6'c < 7r/2. We obtain, 

(2.48) 



LN{x,y) = AN{x,y)- 



{2nf 

and this is the formula we will use. Let 

g{e)^g{d-x)^\og\GN{r,e'';x)\ 



GN{rce'';x)GN{rce'^;y){e'^ 



= -iV log + I log 1 1 + r,e'^ 1 2 + log 1 1 - r,e^« | ^ . 



Taking the derivative gives 
(2.49) 



5'(0)-^sm( 



cos 6 



1-f rce'^pll-rcc*' 



From this formula we see that 

(2.50) 9{e)<g{e,), -n<0<Tr. 
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Taking the absolute values in ( |2.4^ ) gives 



2r, 



(2.51) 



Fix S > (small), li S < 6c < n/2 we can make a quadratic approximation around 
9c to obtain the estimate 



(2.52) 



Jo ~ 



for some constant C. If Be is close to we have to be more careful. A computation 
using ( p.4£ ) shows that there is a constant a > such that 

(2.53) / eS^'^^-a'^'^-Ue < / e-"-^("'-"')'du. 

Jo Jo 

Here we have taken u ~ sin0, Uc — sin 0c- Now, 

/■OO 

(2.54) 



and hence 
(2.55) 



du<Cmm(^^^,^^/J 



^^g{.6;^)-a(6.{i);x) < Cmin( 



sin6'c(e)%/if ' 



1/4)- 



To prove ( 2.54 ) assume that Uc > Set u = Uc(l + s). The left hand side 

becomes 

C 



aKu^s" 



ds = 



On the other hand if Uc < l/K^^^, then the left hand side equals 

cV2 



00 



and we have proved (2.54). Note that sin 6*^(0 = -rT VK^ - - - 1/2)^ 

so ( p. 55 ) can be written 



(2.56) 



es(«'^)-9(«'=(«);^)(i6'< Cmin( 



K^\t{l~t)^{^~l/2)^\' K^' 



1/4 



If ^ > 1/2 + ^J t{l — t) we can take 9c ^ and then ( 2.56|) s till holds by a similar 
argument using ( ]2.49D . Thus, for any Q < x,y < K ,hy ( 2.51 ) and ( 2.56 ), 

|Lw(x, y)\ < CA%{x, y) min(|<(l - t) - (^ - l/2)T'^^ K^") 
(2.57) X min(|t(l - t) - (77 - l/2)2|-i/2^ if i/4)^ 



where 61^(0 is given by (|2.47|) if |^ - 1/2| < ^f(l-t), 61^(0 = if |^ - 1/2| > 
^Jt{l - t) and 



A*N{x,y) = -rcv4Ar(a:,y)e^(''=(«)^-)+s(^=('')^^). 
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Let 5 > and consider x,y such that |^ — 1/2| < y/t{l — t) ~ S and \ri — 1/2| < 
•\/t{l — t) — 5, i.e. we are inside the support of the equihbrium measure. Now, 



= -{z+) = -d+iO 



and a straightforward computation shows that 3fi(i_|_ (^) > 0. Furthermore, if (i_ (^) = 

— -^f{fce'^^)\0=-ea, then d^{£_) — d+{£,). Also, for the ^ and 77 we are considering 
we have that \d+{£,) — d+{ri)\ < C|f — T^j. A standard local saddle-point argument 
now gives, using ^A§i ), ( |2!49| ) and ( |2^ ), 



rcANix,y) 
LN{x,y) = (U 



2tt 



where (a, b — ±), 



0{- 



Write, 
(2.58) 

which defines p in lemma 2.8. Since | eyip{i9c{r])) — exp(z6'c(^))| < G|^ — 77I, 
rfTII) and \d+{0 - d+iv)\ < C\C -ri\,we get 



LN{x,y) 



A^ix^y) 



2TT\d+m 



) <CA*^{x,y){— + \^-7j\) 



K 



which can be written 
(2.59) 

LN{x,y) - A*pf{x:y) simrK{p{^) - p{t])) 

We now investigate A*pf{x,y) and start with the case when both ^ and rj are 
inside the support of the equilibrium measure. Inserting the formulas for k„ and 
w{x) in ( 2.46| ) we obtain 

1 



-K 



1/2 1/2 

y 



(2.60) Aw(a;,y) - (i^-iV+l)| 

Furthermore, a computation shows that 

es(ee(?);-) = r,-^(l + + 2r, cos 0,)"/2(^ ^ ^2 _ 2r, cos ^e)^"^""'/' 
(1 - i)^/22^'/2 



x'^^^iK-x) 



(2-61) iiV/2(l_i)K/2^X/2 

Stirling's formula gives the asymptotic formulas 



{K-x)/2 



(2.62) 

and 
(2.63) 



1 1 

{K - x)K-^x^VK V2^^^(T^1y^^ 



(1 - t)K-^t^VK V27rt(l-t)^^ ^ ^^X^^' 
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Combining ( |I60| ) - (|^) gives 
(2.64) 

A computation shows that 
(2.65) 



A*^{x, y) = /)i^ (l + 0(1) + 0(1^ - ,y|)). 



sin Or 



and hence (2.27) in the proposition follows by combining ( 2.59] ), (2.64) and (2.65). 

We also have to estimate A^(x,y) when ^ or 77 is outside the support of the 
equilibrium measure. Assume that ^ is inside and 77 outside the support (or close 
to the edge of the support), so that 9c{r]) — 0. We have that 

Set 

y/2 / , \ {K-y)/2 



B{y) 



-N 



1 



2r, 



77-1/2 



1 



2rr 



V ~ 1/2 



which corresponds to the expression (|2.6l|) with x = y. We want to show that 

p9(0;y) 

(2.66) 



Biy) 



<7^<1, 



with 7 < 1 if 7/ > 1/2 + \/t{l — t). A computation gives 



es(0;y) 



V/-2 



1 — a 



{K-V)I1 



^7?log^i+(l-,,) log(i5^) 



K 



where a = 1/2 + ^t(l - t), 0<a<77<l. Set 7 = e'''°s f +(i-'7) iog(T-rf ) and note 
that 77 log ^ + (1 — ^7) log(Y5^) < log(a + 1 — a) = by convexity; if 77 > a we get 
a strict inequality. When y is close to K we have to be somewhat more careful in 
estimating the binomial coefficients. We use 



K 



K 



[K - y)^-yyy 



IK 1 1 

— -(l + 0(-) + 0(— - 



:))■ 



Using this formula and proceeding as before we obtain 

t{\-t) ^1/4 



(2.67) 



A*N{x,y) 



-(i + o(. 



instead of ( 2.64 ). Combining ( 2.57 ) and ( 2.67 ) we obtain ( ^.28 ) and the proposition 
is proved. Note that when t = 1/2 some modifications in the arguments above are 
needed. We will omit the details. □ 



2.4. The corner growth model. We can draw the type I DR-paths in a dif- 
ferent way so that they look like the heights of a cascade of discrete polynu- 
clear growth (PNG) models. If we place a W-domino so that it has corners at 
(0, 0), (1, 0), (1, 2) and (0,2) we draw a path by connecting the points (0,1/2), 
(1/2,1/2), (1/2,3/2) and (1,3/2) with straight line segments. We draw a path 
on an E-domino analogously, it is the mirror image of the W-domino in the middle 
vertical line. In this way, from the DR-paths, we obtain height curves, which we 
can think of as graphs of functions hk{x,n), I < k < n, where the A: :th curve 
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goes from Ek — (— n — 1 + fc, 1/2 — fc) to Ak = (n + 1 — fc. 1/2 — fc) in the orig- 
inal coordinate system. We think of hk{x,n) as the height of the level- A: growth 
process at the point x at time t = n. The level-Zc height curve does not intersect 
the level-(fc -I- 1) height curve. Note that the vertical steps always have step size 
±1; hk{x+,n) — hk{x—,n) £ { — 1,0,1}, and the jumps can occur only at points 
X g {-n -l/2 + k,-n + l/2 + k,...,n+l/2- k}. 

A random tiling of the Aztec diamond can be generated by the so called shuffling 
algorithm which we will now describe briefly, see [ p^ and | |3^ for more details. Start 
with Ai. We can tile it by either two horizontal dominoes, with probability 1 — g, or 
two vertical dominoes, with probability q, (compare ( 2.19] ). Assume that we have a 



random tiling of Ak for some fc > 1. We will define a random tiling of A^+i, given 
a tiling of Ak- Call two dominoes which share a side of length two a pair. Two 
horizontal dominoes form a bad pair if the lower one is N and the upper one is S. 
Similarly, a pair of vertical dominoes is bad if the left one is E and the right one is 
W. All other pairs are good. In the first step we remove all bad pairs of dominoes. 
In the second step, for each domino in a good pair we move one unit step, upwards 
if it is N, downwards if it is S, to the left if it is W and to the right if it is E. After 
these two steps what remains to completely fill Ak+i are 2 x 2 -blocks. In the third 
and last step we fill each 2x2 -block with a vertical pair with probability q and 
with a horizontal pair with probability 1 — This generates a random tiling of the 
Aztec diamond Ak+i with the probability ( |2.l[ ), where q = ^^(1 -I- 

The shuffling algorithm translates into a PNG-type growth procedure for the 
cascade of height functions {hk{x,n)}i<ck<n defined above. We will now define the 
discrete PNG-type growth model which we obtain. We have shifted the picture 
1/2 unit upwards compared to the one we obtained from the modified DR-paths 
above. The "height paths" are built from plus steps, which go from (m, n) to 
(m -|- l,n + 1) and consists of straight line segments between the points {m,n), 
(to+ 1/2, n), (to-K 1/2, n-|- 1) and (m + 1, ri + 1), minus steps, which go from (m, n) 
to (to + 1, 71 — 1) and consists of straight line segments between the points (m, n), 
(to + 1/2, n), (m -|- 1/2, n— 1) and (to + l,n— 1). Finally, we have zero steps, which 
are line segments from (m,n) to (to + 2,n). The initial configuration, t ^ 0, has 
zero steps between (n, 0) and (n -I- 2, 0), rt e 2Z + 1. At time i = fc, fc > we do 
the following: 

(i) remove all zero steps; 

(ii) move all plus steps one unit to the left and all minus steps one unit to the 
right; 

(iii) if a plus step and a minus step pass each other in step (ii) they are removed; 

(iv) add zero steps so that we obtain a connected curve from (— oo, 0) to {oo, 0); 

(v) replace each zero step between — (fc -I- 1) and fc -t- 1 with a combined plus and 
minus step independently with probability q. 

We can define a cascade of height curves as follows. The level-TO curve initially 
has just zero steps between (n, —(to — 1)) and (n + 2, —(to — 1)), n e 2Z + 1. At 
each time step we apply the discrete PNG growth procedure independently to the 
levels 1,2, .. . with the condition that the level-TO curve cannot touch or intersect 
the level-(TO -I- 1) curve, to > 1. If that happens in the random growth step, then 
this growth event is suppressed. Note har only a finite number of levels are changed 
at time fc. The shuffling procedure induces an evolution of the modified DR-paths, 
and this is exactly the cascade of PNG growth models just defined. A plus step 
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corresponds to a W-domino, a minus step to an E-domino and a zero step to an 
S-domino. That the level-1 modified DR-paths evolves exactly according to the 
PNG-growth rules is immediate by comparison. All dominoes above the level-1 
DR-path are N-dominoes and move upwards one step. Removing all bad vertical 
pairs corresponds to the annihilation (step (iii)) in the PNG-growth rule. Filling 
in 2 X 2 -blocks is exactly the random growth, step (v) in the PNG growth rule. 
A somewhat more elaborate argument shows that the whole shuffling algorithm 
corresponds to the cascade defined above. We will not give the details. 

Consider the level-1 DR-paths of type I, i.e. the upmost one. In CS-I it goes from 
(1,0) to {n + l,n) through the points {ik,jk), ^ ^ k < p, where (io,jo) — (IjO), 
(ip-jp) = in+l,n) and {ik+i, jk+i) - (ikjk) = (1,0), (0, 1) or (1, 1). The dominoes 
in the north polar zone are the dominoes above this path. Set 

(2.68) n - = max{jfe ; ifc = £}, l<e<n + l. 

Then A = (Ai, . . . , Xn+i) is a partition and this is the partition associated with the 
north polar zone in . (We can see this partition by marking each domino in the 
north polar zone with a point at the center. These points will lie in the integer 
lattice Ci in CS-I.) Set 

A(n) = {(^,J)eZ^; 1 < J < A„ 1 < i < n + 1}, 

which is a random subset of Z^. 

Let w(i, j))(jj)g22^ be independent geometric random variables with parameter 

q, P[w{i,j) = fc] = (1 — q)q'^, fc > 0, and define 

(2.69) G(M,iV)= max V w{i,j), 

IT ^ ^ 

where the maximum is over all up/right paths from (1, 1) to (M,N), see Set 

(2.70) n{n) = {(i , j) G Z^ ; G(M, N)+M + N ~l<n}. 

We get n{n + 1) from Q{n) by independently adding a point with probability 
p = 1 — q to every corner in 17 (n), see |]37| , and because of this we call it the corner 
grow th model. It is proved in [|36[ that if the probability on T(A„) is defined by 
( ^!l| ) with q = w'^{l + w^)^^, then the random sets A(n) and Q{n) have the same 
distribution. 

Let h{x, n) be the level-1 height function in the PNG growth model defined above. 
We can relate the distribution function for this height variable to the distribution 
function of G{M,N). Let Pg = ^(^ + 2, n) be the midpoint in CS-I of the line 
segment from (1,0) to {n + l,n) and set P^, = Pg - fc(l/2, 1/2), fc = 0, ±1, . . . ,±n. 
Assume that both fc and n are even. Then h(k, n) < 2m — 1 if and only if Pk + 
to(— 1, 1) is a point above the level-1 height curve, which happens if and only if 

,n ^ k n + k , „ , x 

{— m + 1,— m+1) e n{n). 

Consequently, for fc and n even, 

_L. 

(2.71) P[/i(fc, n) < 2m - 1] = P[G{— m + 1, — m 1) < 2m - 1]. 
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Using this relation and the asymptotic resuhs for G{M, N) in |^ we can show that 
the fluctuations of the height (and hence of the temperate zone since the height de- 
scribes the boundary of the temperate zone) are of order n^^^ and the ap propriately 
rescaled fluctuations converges to the Tracy- Widom distribution ( 1.11 ). 

As we have seen above the DR-paths can also be related to the zig-zag particle 
configurations. Using this we can relate the distribution function for G{M, N) to 
the distribution of the rightmost particle in the Krawtchouk ensemble. 

Lemma 2.9. If K = t + N + M - I, then 

P\G(M, N) <t]^ Pkt MKn\ max /i,- < t + M - 11. 

' ' i<i<M ■' 

Proof. Considerthe zig-zag path Zr as defined above. It maps to the zig-zag 
particle configuration (/ii, . . . , h^) with hi < ■ ■ ■ < hj.. The relation to the partition 
A defined by ( 2.68 ) is that Xr — n — hp. Hence, G(r, a:) -I- r -I- a; — 1 < n, i.e. \r>x 
in the tiling of An, if and only ii hr <n — x. Thus, by theorem 2.2 , 

P[G(r, x) + r + X — l< n] — Pxr M K q\ max h, < n — x\, 



and this translates into (2.9). □ 

We can now apply the edge scaling result for the Krawtchouk ensemble. Using 
the integral formula for the Krawtchouk polynomials (2.7) and proceeding in the 
same way as in for the Meixncr polynomials, wc can prove that if pt < q{l — t), 
M = [Kt], <t <l, then 



'1/3] 

for each ^ e K.. Here F(^) is given by ( |l.ll| ) and 



(2.72) lim PKr-M.K.,[ max h, < K(3{t) + ^p{t)K'^^] = F(0 



Pit) = (1 - t)p + tq + 2y^pqt{l~t), 

We can now combine (2.9) and (2.72) (allowing a somewhat more complicated 
relation between M,K and t) to give a new proof of theorem 1.2 in |37j. Note 
that in the derivation of ( ^.9|) w e have not used the RSK-correspondence which was 
central to the approach in |37| . 

Let L{a) denote the length of a longest increasing subsequence in a random 
permutation a from Sn where TV is a Poisson random variable with mean a. Then 

P[L{a) <n]= lim P[G(iV, N) < n] 

if we take q = a/N'^, see [|8|. Now, by (p|), 

P[L{a) <n]= lim Pxr .N.n+2N~i,a/N^ [ max hj < n + N - 1] 

A/— >oo ■ ■ l<j<M 
(2.73) = Jim^det(/- K-^^r .N.n+2N-l,a/N^)p{{n+N,...,n+2N-l})- 

Using the formulas (|]^) and (|2.9[) , it follows that the last expression of (|2j|) equals 

det(/ - Pq)£2({„ 

where Ba is the discrete Bessel kernel, 

B {xy)= ^r^ J-('^^)Jy+^i'^^) - Jx+i{2^)Jy{2^) 

x-y 
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and we have rederived a result in Q and | |3q ]. Precise asymptotics for L{a) was 
first studied in Q . We see that the longest increasing subsequence problem can be 
found in a limit of the Aztec diamond. In the same limit the discrete PNG model 
defined above, appropriately rescaled, converges to the PNG model studied in 

2.5. Zig-zag paths for domino tilings of the plane. Consider the squares 
{m,n) + [—1/2,1/2]^, where {m,n) £ 1? . A domino tiling of the plane, which 
we identify with C, is a covering of the whole plane by 2 x 1 or 1 x 2 rectangles, 
dominoes, where each domino covers exactly two of the basic squares. We can also 
think of this as a dimer configuration of the graph with vertices (m, n) and edges 
between nearest neighbour vertices. A domino covers the neighbouring squares with 
centers P and Q if and only if the edge between P and Q are covered by a dimer. 
We will switch between the domino and dimer languages whenever it is convenient. 
Colour the points (m, n) with m + n even black and the other points white, and 
give the corresponding square the same colour. Consider the line y — —x. A 
domino tiling of the plane induces an infinite zig-zag path around black squares 
in complete analogy to the zig-zag paths in the Aztec diamond. We can map the 
zig-zag path to a particle configuration by saying that we have a particle at x if 
and only if the zig-zag path goes from x — 1/2 — i{x — 1/2) to x + 1/2 — i{x + 1/2) 
via a; + 1/2 -I- i{—x +1/2), i.e. an east-south step. Note thst we have a particle at x 
if and only if either the edge from x — \ ~ ix Xo x — ix ox the edge from x — i(x^ 1) 
to X — ix is covered by a dimer. In this way we get a particle configuration in Z. 

There is a unique translation invariant measure of maximal entropy, the 
Burton- Pemantle measure, on the space of domino tilings of the plane, see |^ and 
|45| . This measure induces a probability measure on zig-zag paths and hence on 
particle configurations in Z; we get a point process on Z. We want to show that this 
is a determinantal point process, given by the discrete sine kernel. Let i? be a 
set of disjoint edges, i.e. they do not share a vertex, in the Z^-graph and let Ue i>e 
the set of dimer configurations which contain E. Let P be a white vertex and give 
the edge between P and P + z the weight z, where z — ±l,±i. Assume that the 
edges in E cover the black vertices bi, . . . ,bk and the white vertices wi, . . . ,Wk- It 
is proved in |45j,[Q, using techniques by Kasteleyn, Q], that 

KUe) = aE det(P(6, - u;j))fj^i, 
where is the product of the weights of the edges in E and 

Att^ 2i sin 6* + 2 sin < 

Using this we can prove the following result. 

Theorem 2.10. The probability of having particles at positions Xi, . . . ,Xm in the 

zig-zag point process defined above is 

'sm^{xj - Xk) 



P(x -\- iy) 



(2.74) P\x] = det , , , , 

TT{xj-Xk) y^.fc^i 

Proof. If we have a particle at x, then one of the edges x — ix, x — ix — 1 
or X — ix, X — ix — i is covered by a dimer. We take E^ to be the set of edges 
Xj — ixj, Xj — ixj — Zj, where Zj — 1 or — i, 1 < j < m. Then, qe = zi ■ ■ ■ Zm, 
bj = Xj — ixj and Wk — Xk — ixk ^ Zk, so that 

P{bj - Wk) = P{xj - Xk - i{xj - Xk) + Zk). 
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Thus, 

= X! l^i^E, ) = ^ zi...Zm det{P{xj - Xk- i{xj - Xk) + Zk 

Zj = l or i 2^j=l or ^ 

m 

= ('^) X! n - - - ^<^(i)) + ^<^a)) 

cr&Sm or i j=l 

m 

(2.75) = (a) [] ^(^^i - = det(i^(a;,- - arfe)),%i, 

where 

K{u) = Y zP{u- iu + z) = P{u+1 - iu) +iP{u+ {-u+l)i). 

Zj — 1 or i 

A computation shows that P{—y — ix) = l)^P(a; + iy) and thus 

P(u + 1 - iu) = i(-l)-"-ip(w - z(w + 1)) = -z(-l)"P(w - i{u + 1)). 
We obtain 

K{u) = i{P{u + i{-u + 1)) - (-l)"P(u - iiu + 1)) 



^tt"^ J -TT J -TT 2ism6 + sine/) 



4Tr^J_^J_^ 2isin9 + sin(j> 



-ded(p 



-dedcj) 



Set 



so that 



If we write 



47r2 J_^ J_^ 2i sm{9 - 0) + 2 sin < 



^ ' ' 2tt J_„ 2i sm{e -(t)) + 2sm(f> ^ 



K{u) = -- / e*^"G(6i, u)de 



e^^ +i- (e-'^ + i)z'^ z 



dz 

if. 

we can use residue calculus to see that 
(2.76) G(M) = I^ if-.<^<0 



if < 6' < TT 



Thus, 



^ ' 2ttJ_^ e~'^ + i 2tt 



TZU 



2-KU TTU 
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Inserting this in (2.75|) gives 



p[x] = det ('!!!it(^£_^(_i)^r . det ^'^"5^"^ ~ "'^^ 



TT{Xj-Xk) J^j^^^ V T^{Xj-Xk) /^-fc^i 

□ 

If we consider the zig-zag path through the center of the Aztec diamond, the 
measure on the zig-zag particle configurations has determinantal correlation func- 



tions, (p. 81), by theorem 2.2. Take r = n/2 and p = q = 1/2. In this case the 



equilib rium measure, ( 2.26 ), has density p'{x) — 1/2, < a; < 1. Hence, we can use 



lemma 2.8 to show that the limiting point pricess, as n — *■ oo, has determinantal 



correlation functions with kernel 

K{x,y) = 



sin I (a; - y) 



x-y 



i.e. exactly the same as in theorem 2.10. This is consistent with the conjecture 



|36[ , [|10[ that in the center of the Aztec diamond a random tiling looks like a tiling 
of the plane under the Burton-Pemantle measure. 

3. The Schur measure and non-intersecting paths 
In section 2 we obtained the distribution function for the last-passage random 



variable G(M, A''), (2.69), using the non-intersecting paths in the Aztec diamond. It 
is natural to inquire whether the Meixner ensemble which is used to study G{M, N) 
in p7[ | can also be obtained in a natural way using non-intersecting paths. The 
picture will again be a cascade of PNG-type growth models but different from the 
one studied in sect. 2.4. 

We will define a certain random growth model, or rather a cascade of growth 
models, such that the probability distribution of the heights above the origin is 
given by the Schur measure introduced by Okounkov, The cascade of growth 
models is actually equivalent with the Robinson-Schensted-Knuth correspondence. 
Viennot, [Vi], gave a geometric construction of the RSK-correspondence for per- 
mutations, often called the "shadow" construction, see also |Q, sect. 3.8. A 
permutation in Sn can be described by putting n points randomly in the unit 
square, Hammersley's picture, [0. If one applies the first step in Viennot's con- 
struction and interprets the paths (shadow lines) as space-time paths one gets the 
polynuclear growth model (PNG) as introduced in by Prahofer and Spohn. 
Thus we can equivalently think in terms of a growth model. If one takes the full 
Viennot construction one gets a cascade of growth models as proved by Okounkov, 
|56| . Okounkov did not base his presentation on Viennot's construction, but instead 
used the formulation of the RSK-correspondence in [|j . There is a generalization of 
Viennot's construction to the case of an integer matrix (generalized permutation) 



called the "matrix ball construction", see |26 . By the same argument, this can be 
translated into a growth model and will lead to the Schur measure. This growth 
model is also given in |5^. We will present a somewhat modified version of this 
growth model, which avoids the limiting procedure in The interesting thing 

is that this version leads directly to families of non- intersecting paths, namely the 
standard ones which can be used to obtain the Schur polynomials, [0, Q. 

Let W — {w{i, j)^ be an 71 X n-matrix with non-negative integer elements, 
and set w{i, j) = if n ^ {1, . . . , n}^. The integer- valued height functions in the 
cascade of growth models are denoted by hk{x,t), 1 < k < n, where hk{x,t) is 
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the height above a; e M at time t e N of the level k growth process. The height 
curves x hk{x,t) do not intersect, hk{x,t) — hk+i{x,t) > 1, 1 < fc < n for all x 
and t. The height curves will grow by the addition of unit squares and the growth 
procedure is defined as follows. 

The vertical "sides" of hk{x,t) will be labelled. We can think of the curve x 
hk{x, t) as a lattice path starting at (— fc+1, — 2n+l/2), ending at (— 2n— 1/2) 
and taking unit steps up, to the right or down. Each unit step up is labelled by 
aj for some j, 1 < j < n and each unit step down is labelled by bk for some fc, 
1 < k < n. Call these vertical sides left and right vertical sides respectively. At 
time 0, hk{x, 0) = — (fc — 1), 1 < fc < n and there are no vertical sides. 

Assume that hk{x,t) has been defined for t < rn, — 1, 1 < k < n with labels 
on the vertical sides and such that the distance between a left vertical and a right 
vertical side is always odd. Furthermore hk{x,t) — hk+i{x,t) > 1. We will define 
hk{x, to), I < k < n so that it has the same properties. For each a; G Z, n G N and 
1 < fc < n, we have a set B'''^~^\x,t) of unit squares with labelled vertical sides. 
B^°\x, t) contains w{i,j) unit squares with the left vertical side labelled Oj and the 
right vertical side labelled bj, where i = (t + x + l)/2, j = {t — x + l)/2. Recall 
that w{i,j) = if n ^ {!,..., n}^ so that B'^^\x, t) is empty for odd x at odd times 
t, and for even x at even times t. The B^''~^\x,t), k > 2 are defined recursively 
in the growth procedure. Assume that B'^^^^\x,m) has been defined for some £, 
1 < £ < n. We will define he{x, m) and B^^\x, m). 

(Horizontal growth). Move each left vertical side in he{x,m — 1) one unit to 
the left, and each right vertical side one unit to the right, fill in with horizontal 
segments so that we get a connected curve, and denote the resulting height function 
by /i|(a;, TO — 1). The labels move together with the vertical sides. Since the distance 
between a left and a right vertical side is odd they cannot meet at the same point. 
Set u = he{x — 1, to — 1) — hi{x,m — 1) and v = he{x + 1,to — 1) — hi{x,m — 1). 
If = min(M, v) > 0, then a right vertical side, with labels b^, ■ ■ ■ , br^ (ordered in 
the upwards direction), will cross a left vertical side, with labels a^^, • • -g^. If this 
happens we put z unit squares in B^^\x,m) with vertical sides labelled arj,bsj, 
1 ^ j < -2- This defines B'^^^x, to). 

(Vertical growth). Next we put the labelled squares in ^^^^^^(a;, to,) on top of 
hg{x,m— 1) at X for all x, i.e. between x— 1/2 and x+ 1/2. The result is h^{x,m). 
Note that all vertical sides are labelled, and that by the way that the B^^^^\x, m):s 
were defined, the distance between left and right vertical sides is always odd. 

In this way we recursively define hk{x,t), with labelled vertical sides for 1 < 
k,t<n, X eR. 

Next we want to describe the final configuration. Let A — (Ai,...,A„) be a 
partition and let 7^(A;a) denote the set of all non-intersecting, labelled, up/right 
lattice paths F = {Ffc}^^-^, where F^ starts at (— fc + l,—2n + 1/2) and ends at 
(—1/2, Xj — j + 1) and where all up-steps have an x-coordinate of the form 2{j — 
n) — 1/2 for some j, 1 < j < n. Each imit length of the vertical sides with the 
x-coordinate 2{j — n) — 1/2 is labelled by aj, < j < n. Let V{X;b) be the paths 
we obtain by reflecting the paths in ^(A; a) in the y-axis and putting the label bk 
on each unit vertical side with x-coordinate 2{n — k) + 1/2, 1 < k < n. If we join 
a family of non-intersecting paths from V{X; a) to a family from P(A; b) by adding 
horizontal segments from (— 1/2, Aj — j + 1) to (1/2, Aj — j + 1), 1 < j < n, we 
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obtain a family of non-intersecting, labelled height curves. Let W(A; a, b) denote all 
the families of non-intersecting height curves obtained in this way. 

We claim that H = {hk{x, 2n— 1)}^=! constructed as above belongs to 7i(A; a, b) 
if we put \j = hjiO, 2n-l)+j - 1. We write H = (r, f), where T e V{X: a) and 
r g P{X; b). To see this note that a square with labels aj, bk is introduced at level 1 
at position j—fc at time j+fc—l. In the remaining time 2n—l—(j-|-fc—l) = 2n—j—k, 
the left vertical side moves 2n — j — k steps to the left and the right vertical side 
2n — j — k steps to the right. Thus the Oj-label ends up at a position with x- 
coordinate j — k — 1/2 — (2n — j — k) = 2{j — n) — 1/2, and the 6fc-label ends up at 
a position with a;-coordinate j — k + l/2 + 2n — j — k = 2{n — k) + 1/2. Thus all left 
(right) vertical sides end up to the left (right) of the origin at the correct positions. 
Thus, we obtain a map from the set of n x n integer matrices to [JxH{X; a, b), where 
the union is over all partitions A = (Ai, . . . , A„). This map is invertible, and hence 
we obtain a bijection. That the map is invertible follows from the fact that the 
growth procedure can be reversed! We start at the bottom level, hn{x,t), but now 
right vertical sides move to the left and left vertical sides to the right. If they cross 
we move squares up to the next level, just as we moved them down before. At 
the upper level they tell us where we should split and introduce new left and right 
vertical sides. Those squares that are taken out at the top level give us the entries 
in the matrix. If we take out m sqiiarcs with labels Oj, bk we put the number m at 
position (j, k) in the matrix. (We can also find the position from the time t when 
and the place x where the squares are removed, x = j ~ k, t — j + k — 1.) We 
will not describe all the details of this reverse procedure. In a sense the cascade 
of growth models records the history of the growth process. When two vertical 
segments pass each other at a certain level information is lost at this level, but this 
information is recorded at the next level. 

If the w{j,k):s are independent geometric random variables with P[w{j,k) = 
m] = (1 — ajbk){ajbk)"^, m > 0, then the probability of a particular integer matix 
W={w{j,k))li^, is 

n n n 

n (1 - ajbk) n i^jhr^''"^ = n (1 - ^j^w), 

j,k=l j,k—l j,k=l 



where 



j=i k=i 



If we interpret the labels on the vertical sides of ^(A; a) as weights and define the 
weight, Lj{T), of an element F € P(A; a) as the product of the weights of all vertical 
sides (horizontal sides have weight 1), then we see that the growth procedure defined 
above transports the weights in the correct way; if W has weight ui{W), and W 
maps to (r,r) in W(A;a, 6), then 

oj{W) = a;(r)a;(f). 



The total weight of ah up/right paths from (-fc + 1, -2n+ 1/2) to (-1/2, A^ - j + 1) 
where all vertical sides have a;-coordinates of the form 2{i — n) — 1/2, 1 < j < n, is 
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given by 



fli •■•On — h-Xj-j+kiai, ■ ■ ■ ,an); 

Tni-\ l-rrLn — Xj 



We have vertical steps with ^-coordinate 2{i — rt) — 1/2 and these have weight 
tti. Here /im(ai, . . . , a„) is the complete symmetric polynomial of degree to in n 
variables; hm{a) = if to < 0. 

The Lindstrom-Gessel-Viennot method, which is the discrete analogue of the 
Karlin-McGregor theorem, see for example p. 98 for a precise statement, gives 

(3.1) ^(r) = det(/iA,-,+fc(a)),Vi =*A(a). 

The expresion in the middle can be taken as the definition of the Schur polynomial 
sx{a). The growth procedure above defines a map S{W) — A from the integer 
matrix W to the partition A defined by the succesive heights. We obtain 

n 

(3.2) P[S{W)^X]^ n (1 - «^-^^)'^(^) 

n 

= n (1 - E '^(r)^(f ) 

j,k=i rev{\;a),rev{\-b) 



(1 - Qjbk 

3,k=l 



Sx{a)s\{b) = Psclmr[A], 



the Schur measure on partitions, introduced in p5[ . 

Remark 3.1. Note that an element F — {TkY^^i ^(-^5 a) corresponds to a 
unique semi-standard tableaux T with shape A, sh (T) — A. If has rj vertical 
steps at 2{j — n) — 1/2, we put Vj j:s, j — 1, . . . ,7i in weakly increasing order 
in the Xk boxes in row k. Similarly an element in V{X]b) gives a semistandard 
tableaux D. Thus we obtain a one-to-one map W {S,T), which is exactly the 
RSK-corresponde nce. If we let rrijiT) denote the number of j:s in T, we see that 
we also have, by ( |3.l|) , 

reV{\;a) T:sh=X 



which is the combinatorial definition of the Schur polynomial, see 



Consider now the random variable G{M, N) defined by ( 2.69 ). From its definition 
it is clear that we can compute G{M, N) recursively by 

(3.3) G(M, N) = max(G'(M - 1, N),G{M, N - 1)) + w{M, N). 
We want to show that, for 1 < M, N < n, 

(3.4) G{M,N) = hi{M-N,M + N-l), 

in particular G(n, n) = /ii(0, 2n — 1) = Ai, which is a well known property of the 
RSK-correspondence, Q , Gsl . It follows from the definition of the growth process 
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above that 
(3.5) 

t + X + 1 t — X + 1 , 

h{x, t) = max(/ii(a; — l,t — l),h{x,t — l),h(x + — 1)) + w{ , ] 

(Recall that w{j, k) ^ if {j, k) <^{1,..., n}^.) Write G{j, k) = h{j -k,j + k- 1). 



Then ( |3.5| ) becomes, 

(3.6) G(j, k) = max(G(j - 1, k),G{j ~\,k~ \), G{j, k - 1)) + w{j, k). 

We will now use induction, li M + N - I = I, then M = N = I and G(l, 1) = 
w(l, 1) = G(l, 1). Assume that G(M, N) = G{M, TV) for M + - 1< k. We want 
to show that it is true for M + — 1 = fc. Note that, by (3.6), 

(3.7) 

G(M - 1, iV - i) = max(G(Af ^\,N ^ \),G{M - 1, iV - 1), G(Af -\,N ^ ^)) 
since w{M — 1/2, N — 1/2) — 0. By our assumption and (|3 

(3.8) G{M - 1, iV) = G(Af - 1, iV) > G(A/ -\) 
and 

(3.9) G{M, N^l) = G{M, N - 1) > G{M ^). 
Furthermore, by our assumption, G(A/ - 1, TV - 1) = G(A/ - 1, iV - 1) and G(Af 



1,A^) > G{M - 1, A^- 1) by (|3J). Combining this with (^, (^ and (^ we 
see that 

(3.10) max(G(Af- l,Ar),G(A//,Ar- 1)) > G{M -^,N -^). 

Consequently, by (3.6), our assumption and ( 3.10| ) 
(3.11) 

G(Af, A^) = max(G(A./ - 1, Af), G(Af -^,N- i), G(Af, N - 1)) + w{M, N) 
= max(G(A/ - 1, N),G{M, N - I)) + w{M, N) = G(Af, A^), 

which co mpl etes th e pr oof. 

From (3^) and (3^) we obtain 

n 

P[G{M,N) < - [ n (1 - sa(«)^a(6). 

j,k=l A:Ai<t 

Using the fact that the Schur measure has detcrminantal correlation functions, 
which was proved in |Q , see also Q , we see that this equals a Fredholm determi- 
nant with a certain kernel, and this can be exploited for the asymptotic analysis. 

Remark 3.2. The case when w(i, j) are independent exponential random variables 
with mean 1 can be obtained as a limit of the geometric case as discussed in . 
We can take the same limit in the construction above and this leads to a continuous 
analogue of the RSK-correspondence. The resulting picture of paths can be viewed 
as two families of n non- intersecting Poisson processes with rate 1/2. (We take 
Oj = 1 — fjlL^ bj — 1 — 1/2L, so that ajbk ~1 — 1/L when L is large. This gives 
the rate 1/2 for the limiting Poisson processes on both sides.) Take the vertical 
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axis in the negative direction as time axis, and the horizontal axis as counting the 
number of events in the Poisson process. If the heights above the origin in the 
cascade are ii > • • • > i„, then the A::th process Ak on the left (and right) start at 
1 at the time —tk and end at n + 1 — fc at time 0. The probability for Xi to go from 
to n — j in a time interval of length ti is 

(3.12) P[X,{U) = n-j] = e-*-/2M)r!. 

(n~ J)] 

The Karlin/McGregor theorem can be generalized to unequal starting times, see 
|42[ , and we find that the probability for the n non- intersecting paths to the left 
with the specified initial and final positions is 

(3.13) det(e-*'/2^^l/^)i:,.^i = 2-"("+i)/2 f[ ^A^it) f[ e'''/\ 
The possible heights lie in [0,cxd), so we obtain the probability density 



1 " 

(3.14) ^A^{t)'l[^' 
where 



P n 



Hence we have rederived the result of proposition 1.4 in (in the case M = N) 
using non-intersecting paths. The probability density ( 3.14[ ) is a special case of 



the Laguerre ensemble. It is also possible to consider the case when w{i,j) is 
exponential with parameter ai + j3j . This leads to a continuous analogue of the 
Schur measure, which can also be obtained as a limit of the Schur measure defined 
above with ai = 1 — Ofi/L, hi = 1 — (3i/L as L — > oo. Using the methods of | |39| we 
can derive the correlation functions for the continuous Schur measure, and we can 
also obtain proposition 1.4 in in the case M ^ N. 

4. Random walks and rhombus tilings of a hexagon 

4.1. Derivation of the Hahn ensemble. In this section we will explore the 
relation between non-intersecting random walk paths and another tiling problem. 
We will consider random tilings of an abc- hexagon with rhombi, see p^ , which 
are directly related to so called boxed plane partitions, Js^ . An abc-hexagon has 
sides a, 6, c, a, b, c (in clockwise order) and equal angles. We want to tile this region 
with unit rhombi (often called lozenges) with angles tt/3 and 27r/3. The number of 
possible tilings is given by MacMahon's formula, |Q, 



(4.1) ^(a,6,c)^nnn 



i=ii=ii=i •< 

We obtain a random tiling by picking each tiling with equal probability. It is well 
known that a rhombus tiling can be described by non-intersecting random walk 
paths, see for example Q and below. This is the approach we will adopt here. We 
will give a new derivation of the Hahn ensemble introduced in . 

Let a, 5, c > 1 be given integers. Take e = (0, 1/2) and f — (-\/3/2,0) as basis 
vectors in our coordinate system; all coordinates will refer to this choice of basis 
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vectors. Our hexagon has corners at Pi = (0, 0), P2 — {b, —b), P3 — {a + b,a ~ b), 
P4 — (a + b, a — b + 2c), P5 — {a,a + 2c) and — (0, 2c). Consider random walks 
S''{m), < m < a + b, starting at (0, 2(fc-l)) and ending at (a + 6, a-6 + 2(fc- 1)), 
1 < A: < c, 

m 

where = ±1 are independent Bernoulh random variables taking each value with 
probability 1/2. Assume that these random walks are non-intersecting. We restrict 
to the case a> b; the case a < 6 is completely analogous. Set 

— m if < TO < & 
m — 2b if b < m < a 
TO — 26 if a < m < a + b, 

TO-|-2(c-l) ifO<m<6 
TO -|- 2(c — 1) if 6 < m < a 

2a — TO + 2(c — 1) if a < m < a + b, 

and 7,„ = (/?,„ - ajn)/'^- Then a,„ < S''{m) < (3rn- Set Xk = {S''{m) - a„i)/2 
and note that < a;^ < 7™, 1 < fc < c. These numbers describe the points where 
the random walks intersect a vertical line. Note the analogy with the previous 
problems. Think of the positions of c particles in a discrete gas confined 

to {0, ... , 7„i}. Let ^1 < • • • < £,L^, 

if < m < 6 
if 6 < m < a 
m if a < m < a + b, 

be the positions of the holes. The holes correspond to the positions (to, a™ + 2^^) 
in our coordinate system. 

There are three types of rhombi. Type I which are spanned by e + f and 2f , 
Type II which are spanned by e — f and 2f and Type III (called vertical) spanned by 
— e -I- f and e -I- f . (We will sometimes call type I and II horizontal.) Given the non- 
intersecting random walk paths we can now tile the hexagon with rhombi as follows. 
If X|!^_^_l = 1 we put a type I rhombus at (to, S''{m)), and if Xm+i = —1 we put a 
type II rhombus at (to, S''{m)). Finally we put a type III rhombus at each of the 
hole positions (to, am+2£,k)- Note that the vertical rhombi are associated with holes 
in the gas, whereas each particle is associated with a horizontal rhombus. Our non- 
intersecting random walks correspond to picking all of the possible hexagon tilings 
of the abc-hexagon with equal probability. We want to compute the probability 
measure induced on the particle/hole configurations on the vertical line x = m. In 
order to be able to formulate the results we first define the Hahn and associated 
Hahn ensembles. 

The Hahn ensemble, js^ is a probability measure on {0, . . . , iV}™ defined by 
(4.3) P]^i\h] = ^A„,(/i)2n^(,-« (/.,), 



13m = 



(4.2) 



7™ + 1 - C 
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where a, f3 > —1 are given parameters, 



(4.4) 

a weight function and 
(4.5) 



"at 



it) = 



{N + a ~ ty.jf] + ty. 
t\{N^ty. 



,(",/3) 



he{0,...,N}'^ 3 = 1 

a normahzation constant. This ensemble is related to the Hahn polynomials, 



which are orthogonal on {0, ...,-/V} with respect to the weight (4.4). For some 
facts about these polynomials see the proof of lemma 4.2 below. Using the leading 



coefficients of the normalized Hahn polynomials a standard computation, |53], gives 



(4.6) 



_ ™, 'fr' + MP + Ma + f3 + j+N + l)!(a + (3 + j)! 



n 

J=0 



(a + /3 + 2i)!(a + /3 + 2j + l)!(iV-j)! 



The associated Hahn ensemble on {0, . . . , iV}™ is defined by 



(4.7) 



where Z 



^ m 

H=-7^A„,(/^)2n 



(",/3) 
N.m 



-(",/3) 
"at 



Na 



is a normalization constant and 



(4.8) 



~(a,/3) 



W 



N 



it) 



t\{N ^ty.{N + a~ty.{(3 + ty 



The ensembles (4.3) and (4.7) are related by a particle/hole transformation, com- 
pare ( p. 17 ) above, and see the proof of theorem 4.1 below. 

Let Pm{xi, . . . ,Xc) denote the probability of having the particle configuration 
xi^ . . . , Xm at time m (along the vertical axis x = m), and P(^i, . . . , ^i^^), 
given by ( |4.2| ) , the probability of having the hole configuration , . . . , £,l^ at time 
m. 

Theorem 4.1. If a,b,c> 1 are given integers, a >b, and we define a„i = |a — m\, 
Kn = \i> — m\, then 



(4.9) 

and 

(4.10) 



^'m(ei,---,a„ 



Proof. The number of random walk paths from j to fc in m steps is (m + k — j 
even) , 



m 

ra+k—j 



e^+.-,(l'"). 



where 1™ = (1, . . . , 1) G N™ and e„(x) is the elementary symmetric function. We 
can now use the Karlin- McGregor, Lindstrom-Gessel-Viennot argument to see that 
the number of non-intersecting random walk paths from (0, 2(k — 1)) to (m, 2xk + 

am) is 

(4.11) AM = det{es^+,,^,+,iin)lk=u 
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where 6m = + CKm)/'2- Introduce the shifted particle positions 

Sfc = Xc+l-k + e {0, . . . , 7m + Sm}, 

which means that we have introduced Sm extra holes at the positions j — 1, 1 < 
J < 5m. Define the partition A by 

— Sk + k — c, 1 < k < c. 



By reversing the order of the rows and columns in (4.11) we obtain 
(4.12) = det(e,,_,.+,(l")),^^,^i = .v(l"), 



where A' is the conjugate partition to A and s\ is the Schur polynomial, see (3.4) 
above and pO[ |. 

If we set rj = c + j - 1 - X'j, 1 < j < Lm + Sm, then {ri, . . . , tl^+s^} U 
{si, . . . , Sc} — {0, . . . , 7„i + Sm} so the give the positions of the (shifted) holes 
including the extra holes. The positions of the original holes are given by 

(4.13) - rs^+j ~ Sm, 1 < J < Lm- 

Note that L.,„ + 5m = rn if < m < a and Lm + Sm — a if a < m < a + b. Let 
fi — A' = (A'^, . . . AJ„) if < m < a, and n — (A'^, . . . A^, 0, . . . , 0) {m — a extra 
zeros) a a < m < a + h. Then, 



(4.14) A™(x)=sa'(1™) = Sm(1™)= n 

1< i<j<7n 



i - i 



by the classical formula for a Schur polynomial. We now want to rewrite the right 
hand side of (4.14) in terms of ^ using (2.17) and ( 4.13| ). Some computation gives 

(4.15) Am{x) - C(m, a, 6, c)Ai„ (0, 
if < m < 6, 

(4.16) Am [x) - C(m, a, 6, c) J] fe'+J^'^^^ A^^ (C), 
if 6 < TO < a, and 

(4.17) Amix) . C(TO, a, 5, c) H (0 + - + c - 1 - e.) 
if a < TO < a + 6, where C{m, a, 6, c) is a constant, e.g. 

(4.18) ^(-'«'^'^)= n"(7T7b^ 

for b < m < a. 

We will also write Wm(C) = ^m(a;) = the number of non- intersecting random 
walks ending with hole configuration ^ at time to. The number of possible non- 
intersecting random walks coming from the right side and going in the other direc- 
tion is Wm' {^'), where = c + Lm — 1 — ^l^+i-j and m' = a + b — m. This follows 
from the symmetry of the hexagon. The total number of tilings, given ^, at time 
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m is then WmiO^m' {£.')■ This can be computed by using ( 4.15 ) - ( 4.17 ) with the 
result 

(4.19) W^iOWm'iO = C*(m,a,5,c)Ai„(0'n<r'''"HCj), 



where C*{m, a, h, c) is a constant. Using ( 4.18| ) we see that, for 6 < m < a 

6-1 ^ 

(4.20) C*(m,a,6,c) = C(m, a, 6, c)C(m', a, 6, c) = TT — — — - 

'~_Q (j + "T- ~ oj!(j + a — to)! 

It follows that the total number of tilings is 



Ce{o,...,7,^}^ 
C*{m,a,b, c)Z^^2lZ' ^ 



by (4^). This proves ( 4.10 ). Note that by combining ( 4.2C ) and ( fl.q ) we obtain, 
after some computation (where we take m = b), 

j!(a + c + j)! 



N{a,b,c) = l[ 



3=0 



(a+j)!(c + j)!' 



which proves MacMahon's formula ([4.11) . 

We now want to go fro m th e variables ^ to x. Since {xi, . 
{0, . . . , 7m} we can use ( 2.17 ) to get 



.(0 



n^l n 



\j=i 



..,Xc}U{^l, 

Ac(x). 



Using this it is straightforward to show that 



WmiOWm'iC) = a(™,a,6,c)Ae(x)2 ]J^i)^"j^>'''")(x,), 



and (O) follows. 



□ 



4.2. Some asymptotic results. Random tilings of a hexagon shows the same 
type of arctic ellipse effect as the Aztec diamond, see .... We will have polar zones 
associated with each corner of the hexagon. Consider the vertex Pq — (0, 2c) of the 
hexagon. We will say that two type I rhombi are adjacent if they share an edge. The 
polar zone associated with Pq is now defined as follows. If there is no type 1 rhombus 
Ro having Pg as a vertex, the polar zone is empty. Otherwise a type 1 rhombus R 
belongs to the polar zone if there is a sequence of rhombi Ro, . . . ,Rk = R such that 
Rj and Rj+i are adjacent. Consider the horizontal rhombi immediately to the left 
of the line x = m, the m:th column, there are c of them. If Zm — Cl„ = max^j 
is the position of the last hole, i.e. vertical rhombus on the line x = m, then all 
the horizontal rhombi above it in the to :th column are of type I and belong to the 
polar zone of Pg. Hence, the boundary of this polar zone is obtained by joining the 
points Ai, Bi, . . . ,Aa, Ba, where Am = {m-l, Um + '^Zm + l), B„i = (to, Q!m + 2Z„) 
with straight line segments. The boundary of this polar zone is thus related to the 
position of the rightmost particle in the Hahn ensemble. 
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The asymptotic position of the rightmost particle and its large deviation prop- 
erties can be investigated using the results of |3^, sect. 2.2. Consider the Hahn 
ensemble (4.3). Assume that m/N ^ i € (0, 1) and -^(a, /?) (aoi Po) b,s N ^ oo. 
We have the limit 

V{s) = - lim -W^"^^\ms)^--U{ts), 
N-*oo m t 

where 

U (s) = (1 + ao - s) log(l + ao- s) + (/3o + s) log(/3o + s) 
— slogs — (1 — s) log(l — s) — C. 
Here we have introduced the modified weight, 

\ N J 

which only differs by a multiplicative constant. The equilibrium measure^ ui*q"°'^"^ (s)ds, 
for the Hahn ensemble is the unique solution of the constrained, weighted varia- 
tional problem 

/ .i/t .i/t .i/t \ 

(4.22) Fy=inMy J \og\a ~ s\-^u{a)u{s)dads + J V{s)u{s)dsj, 

where A = {u e L'^[0,l/t] ; J^^* u = 1 and < it < 1}. From iQ, theorem 2.2, 
we obtain the following large deviation result. Let R — R{t, ao, f3o) be the right 
endpoint of the support of the equilibrium measure, and let e > be given. There 
are functions L{R — e) and J(i? + e) such that 



(4.23) hm i^logP^"„f [Imax/j, < i? - e] = -2X(i?- e). 



m 

and, if J{R + e) > for e > 0, then 

(4.24) lim log P^";f ' [— max h,>R + e] = -2J{R + e). 

We always have L{R — e) > if e > 0, but we must prove that J{R + e) > if 
e > 0. The function J is defined by 

J{x) = inf g{x), 

where 

(4.25) 9i^)^J^ ^og\x -y\~^u,q{y)dy + -V{x) + -J^ V{y)u,^{v)dy ^ Fy . 

By the general theory for the constrained variational problem ( 4.22| ) , see , g{x) > 
for x> R. Now, for x > R, 

9"{x)= r -!^sM^dy+iv"{x) 

Jo [x - yV 2 

and V"{x) = -tU"{tx) with 

"0 



U"{s) 



(l + ao-s)(l-s) (/3o + s)s' 
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SO V"{x) > and g is strictly convex. Consequently g{x) > ii x > R, which is 

what we wanted to prove. 

Note the asymmetry in ( 4.23| ) and (4.24). Just as for the arctic circle large 



inward fluctuations of the temperate zone have much smaller probability than large 
outward fluctuations. 

We will now compute (a part of) the arctic ellipse. For simplicity we restrict to 
the case a — b. The general case can be handled similarly but the computations are 
somewhat more involved. One approach is to compute the equilibrium measure, and 
hence its suport which gives the arctic ellipse, by solving the variational problem 
( 4.22| ) as was done for the Krawtchouk polynomials in |l4j . Below we will instead 



use the approach of |50 which is based on the recursion formula for the Hahn 
polynomials. 

Lemma 4.2. Let R{t,ao) be the right endpoint of the support of u^eq"°'°'°'^ o,s de- 
fined above. Then, 

1 .1 1 



(4.26) R{t, ao) = - sup (- + -— ^ s{l - s){s + 2ao){s + 2ao + 1)). 

t o<s<t 2 2{s + ao) 

Proof. Let q„ ~ (1^'^n\'^) denote the normalized Hahn polynomials which are 



orthonormal on {0, . . . , N} with respect to the weight ( 4.21 ), 



quN \x) = -r^M-n, -.T, n + a + /3 + 1; /? + 1, -iV; 1) 

an,N 

(4 27) _ (-1)" ^ (-n)fc(-x)fc(n + g + /3 + 1)^ 



where dn^^ > and 
dl 



duM ^ {f3 + l)k{-N)kk\ 



[a + 13 + l)(a + l)n[N + a + /3 + 2), 



(^)(2„ + a + /3 + l)(/3+l)„(a + /3 + l). 

The leading coefficients are 

^ (-1)" (n + a + /? + l)n 
(4-2») Kn,N = / _^^\ I at\ ' 

dn,N [a + l)n(-N)n 

The polynomials satisfy the recurrence relation 

Xqn — anMln-l + b^Mln + O-n+LN In+l^ 

where 



_ n(n + a)(n + a + (3 + N + I) [N - n + l){2n + a + (3 + 1){I3 + n){a + (3 + n) 

0-n,N 



{2n + a + [3){2n + a + (3 + l)V [a + n){n + N + a + (3 + l)n{2n + a + [3 + I) 
and 

+ /3 + l)(7i + /? + 1){N - n) n{n + a){n + a + (3 + N+l) 



bn,N — 



N{2n + a + f3+ l){2n + a + P + 2) N{2n + a + (3){2n + a + /3 + 1) ' 



Consider the rescaled equilibrium measure, u%^{s) = iwoq(s/i), < s < 1- The 
paper ||5^ tells us how to compute the support of u*q (and also the measure itself) 
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using the asymptotics of the recursion coefficients. We restrict to the case ao — Pq. 
Then, n/N t e [0, 1), 

lim an,N = a{t) = ^^(1 - t){t + 2ao){t + 2ao + 1) 

N^aa 4(t-|-Q;oj 

lim bn,N = b{t) = i. 

iV— too 2 

According to |M, p. 171, the right endpoint of the support of u*, i.e. tR{ao, t) is 



tR{ao, t) = sup {b{s) + 2a{s)). 

0<s<t 



We can now formulate the result we obtain for the arctic ellipse. 



□ 



Theorem 4.3. Consider the abc-hexagon, assume that a — b and rescale the size 
of the hexagon by a factor l/c. We let the size of the hexagon grow in such a way 
that a/c ^A>0 asc— >oo. Pick the ON- coordinate system in which the limiting 
hexagon has corners at ±(— -^A, i), ±(0, ^(1 + A)) and ±(^A, i). Let X be the 

intersection point of the line x = t, — <r< —-^3^, with the inner part of 
the boundary of the reseated polar zone at Pg. Then, 



2 



(4.29) X \/2A + 1^/1/4 -t2/3A 

in probability as c —i- 00. For any e > there are constants /(e) > and J(e) > 
such that 

1 



(4.30) - logPlX > V2A + lv/l/4-rV3A2 + e] ^ -/(e) 

c 



an I 



d 



1 



(4.31) logP[X < x/2ATTv/1/4 - t73A2 - e] ^ - J(e) 
as c ^ 00. 

Proof. We are in the case when 0<m<a = b, son = L„i = m, N = 7„j = 
m + c^l, a — a — m — P^b^m. Assume that m/c— >/i>0 as c—>oo. 
Then t = /i(l + fi)-^ and ao = (A - ^)(1 + fi)^^ = (1 - i)A - t. Consider the 
points Bm — (ni, am + 2Zm) which describe the inner boundary of the polar zone; 
a-m = — m. We see that 

(4.32) Z,n/Lm^ Rit,ao) 

in probability as c — *■ 00 Here we use the large deviation formulas ( 4.23| ) and ( 4.24 ) 
together with 4.1. To get an ON-system (with the sides of the rhombi = 1 we have 
to rescale the coordinates to /?,„ — (-^m, — ^ + ^m)- Thus, 

c Z c 2c c Lm 2 2 

We also have to translate the coordinate system so that we get the origin at the 
center of the hexagon. We then get the coordinates 

(4.33) (^(^_A),-^+A./?) 



38 



K. JOHANSSON 



for the limit of the point X on the arctic eUipse. Now, 



(4.34) 



R{- 



A 



7i + 1 
for < < A(A + 

1 



1 



1 



1 



2^ 



3('S) = + 



To see this we use lemma [4.2 
1 



h 1)m(2A 
and set 



m) 



Vs(l-s)(.5 + 2ao)(s + 2ao + l) 



2(s + ao) 

with ao = (1 — 1)\ — i. We must have g{s) < 1 because the support is restricted to 
[0, A computation shows that g{s) is strictly increasing in [0, sq] and strictly 
decreasing in [sq, 1] for some sq & (0, 1) For <t < sq, 



sup g{s) = g{t) = - + — — 

0<s<t ^ 



2A 



2A + 1 



(4.35) 

0<s<t ^ V -r J- 

where we have inserted ao = (1 — t)\ — t. Suppose sq < A(2A + 1)^^. Then 
g{so) > g(A(2A + 1)-^) = 1, which is impossible. Thus ( |4.35| ) holds for < t < 
A(2A + l)"i and ( |4.34D follows from ( |4.26| ) since t = /i(/i+l)"^ Set r = ^(^-A) e 

/_V3 _V3 



2 A+1 



-). Then the Imiting point (4^) becomes (r, V2A + 1^1/4 - t2/3A2). 
The large deviation formulas ( |4.30| ) and (4^) follow from ( [f.23| ) and ( [4.24[ ). □ 

The 1-dimensional marginal probability in the Hahn ensemble (E|^) is, fs^. 



(4.36) 



,(",/3) 



^ n— 1 



(",/3)/ 



fc=0 



for t e {0, . . . , iV}, where q["jsf ■* are the normalized Hahn polynomials ( 4.27 ). Hence, 
the probability of finding a hole at position i, i.e. the 1-point correlation function, 
is nu^^'^\t), and consequently the number of rhombus tilings of the abc-hexagon 
with a vertical rhombus at position t on the line a; = to is 



(4.37) 



it). 



k=0 



by theorem 4.1. (This can be rewritten using the Christoffel-Darboux formula.) If 
we use the explicit formul a (|4.27 ) for the Hahn polynomials, we obtain an explicit 
formula for the quantity (4.37). This quantity has been investigated in [|8l, p3. 



25 1 . We will not attempt to rewrite (4.37) in oeder to compare it with existing for- 



mulas. We can also consider the number of tilings with vertical rhombi at specified 
positions ti,...,tr on the line x — m. This is given by N(a,b,c) times a deter- 
niinantal correlation function like ( ^.1C| ) but where we now have instead the Hahn 
kernel given by the formula ( ^.11 ) with the Hahn polynomials and the Hahn weight 
instead. It follows from the general theory in |^ that the 1-dimensional marginal 
probability converges weakly to the equilibrium measure. This should also hold 
pointwise, i.e. 



(4.38) 



lim u 

N^oo 



(a,/3) 
JV,n 



(M) 



if jj{a,(3) (ttOj/^o) SiJid n/N t a.s N ^ oo, but we do not have good enough 
control over the asymptotics of the Hahn polynomials to prove it. The correspond- 
ing result for the Krawtchouk ensemble follows by the same methods as was used 



to prove lemma 2.8 
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As mentioned in the beginning of this section a rhombus tihng of an abc-hexagon 
can also be interpreted as a boxed planar partition where the sides of the box are 
a, b and c. We want to relate the height of the planar partition surface above a 
certain point to our particle configurations. Let {x,y,z) be the coordinates in the 
planar partition coordinate system (the cc-axis through Pi, the y-axis through P3 
and the z-axis through P5 after projection). Consider the line x = m. We start 
from the point (r^, Sm, 0), where 

if < m < 6 
m ifb<m<a + b 



and 




I m if < m < b 
lo nb<ni<a + b 

As we go along the line a particle means that the z-coordinate is increased by 1, 
whereas the x- and y-coordinates are fixed. A hole means that the z-coordinate 
is fixed, but the x- and y-coordinates are reduced by 1. Let Xm{k) denote the 
position of hole number k. This hole corresponds to the position (r^ — fc, Sm — k) 
in the xy-plane and the surface height above this point is equal to the number of of 
particles in {0, ... , Xm{k)}, which equals X,n{k) — fc + 1. Thus, the planar partition 
height function is given by 

(4.39) H{r^-k,s^-k)^X^{k)-k + l. 

Let Ym{n) denote the number of holes in [0,n\ in the particle system on x = to. 
Then 

(4.40) P[X,nik) <n] = l- P[y,„(n) < k]. 

Assume that m/c —ffi>0, L„i/j„i — > t > 0, n/L„i r and j^^{\a — m\, |5— to|) — > 



(ao, Po): as c — > cxd. Then theorem 4.1 and the general theory of discrete Coulomb 
gases in [B^ shows that 



1 



L 



(4.41) ^E[Y„,{n)]^ I ui*,"<"^"Hs)ds 







as c — > 00. Also, it is possible to prove large deviation results, analogous to those in 
Q for the semi-circle law, in this case too. Using this and ( 4.39| ), (4.40) and (4.41) 



it is possible to compue the asymptotic shape of the planar partition surface (law 
of large numbers) as well as large deviation results. See jl^ for this type of results 
proved by other methods. The asymptotic shape can be computed if we know the 
equilibrium measure. 

From ( [1.39D and (|!4^) it is clear that the surface fluctuations are directly related 



to the fluctuations of the number of paticles in an interval in the Hahn ensemble. 
Hence, in analogy with the Krawtchouk case, the surface fluctuations should be 
Gaussian with variance proportional to logc as c — > 00. This could be proved 
provided we had the same control of the Hahn kernel as we have of the Krawtchouk 



kernel in lemma 2.8, but this remains to be done. 

We will end this section by showing one way of finding CUE in a random rhombus 
tiling of an abc-hexagon in the limit as the size of the hexagon goes to infinity. 
Consider the m holes in the TO:th column, 1 < m < b = a. we want to compute 
the probability distribution of these to holes, with to fixed, as the hexagon grows. 
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Let £,1, ■ ■ ■ ,£,m be the positions of the holes. By theorem 4.1 the probabiHty of this 
configuration is 



(4.42) 



^ m 
la-ni.a^ni) A™('?)^ II 



(7m + a - m - £j)l{a - m + £j)\ 



6!(7m-6)! 



where 7,„ = to + c — 1. Denote the corresponding expectation by E'^^^^m 



m.a—m 1 



Proposition 4.4. Let c — > oo, a/c — > A > and keep to > 1 /iied. Let / 
he a continuous, bounded, symmetric function. Then, 

lim 4^;™.-™)[/(il^^!^, . . . , ^"-'Z"'/^ )] 



(4.43) 



1 /■ " 2A 2 

—-- / A„(a;)2lle"^'=?/(xi,...,x,„)d™x, 



where Zm{X) is a normalization constant such that the right hand side is 1 when 

Thus, in the Hmit, the positions of the holes (vertical rhombi) on the TOith vertical 
column are described by to x GUE. 



Proof. The proof is straightforward. By (4.42) the expectation in the left hand 
side of ( 4.43 ) equals 

1 



j{a—m,a—rn) 



E 



?e{o,...,7™}" 



Am(OV(^^,...,^^) 



i=i 



(7 + a - TO, — ) ! (a — TO + ) ! 



In this expression we write — nj -\-^ /2 and use Stirling's formula to approximate 
the factorials. This leads to an expression which is a Riemann sum. The normaliza- 
tion constant is the same Riemann sum but with / = 1. After cancelling common 
factors we see that the remaining quotient of Riemann sums converges to the right 
hand side of ([f.43|). □ 



5. A DIMER MODEL ON A CYLINDRICAL BRICK LATTICE 

In this final section we will consider a dimer model on a hexagonal graph on 
a finite cylinder. The graph is sometimes referred to as the brick lattice, |6^. A 
dimer covering of this lattice can also be thought of as a certain cylindrical rhombus 
tiling, see p5| . The dimer covering has an equivalent description in terms of non- 
intersecting random walk paths, [p^ . In contrast to the previous non-intersecting 
path models we do not have a fixed number of paths. We will again be interested 
in the point process we obtain by looking at where the random walks are at a fixed 
time, and the analysis is based on the methods of 

The graph Gm,n which we will consider is defined as follows. The vertices are 
Vj^k = (-1/2 + j,k), < j < 2M - 1, < k < 2N and we obtain a graph on 
a cylinder by identifying Vjj^ and Vj^2m,k 

for all j, k. We have vertical edges be- 
tween Vj,k and Vj^k+i, and horizontal edges between V2j.2k and V2j+i,2k and between 
V2j+i,2k+i and V2j+2,2k+i- A dimer covering of Gm^n can equivalently be described 
by non-intersecting random walk paths, |l^,(53|. In the same coordinate system 
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as was used to define Gm,n we have simple random walk paths S{t) with steps 
±1, with initial position S'(O) € {2A: ; < /c < A^} and which arc required to sat- 
isfy < S{t) < 2N. Since we have a graph on a cylinder we and the paths to 
live on the cylinder also so we require that S{t + 2M) = S{t). We have L non- 
intcrsccting paths, where < i < A^. (Actually the condition < S{t) < 2N 
implies that L < iV — 1.) Let Pm,n be th set of all such families of non-intersecting 
paths. We now describe the 1-1-correspondence between Vm,n and the set of all 
dimcr configurations on Gm.n- A dimcr covers a vertical edge if and only if wc 
have a random walk step, in one of the paths, which intersects the vertical edge. 
A dimer covers a horizontal edge if and only if no random walk path hits this 
edge. If the horizontal edge V2j.2kV2j+i,2k is not covered by a dimcr, then one of 
ei = V2j+i,2kV2j+i,2k+i 01 62 = V2j+i,2kV2j+i,2k-i must be covcred. A random walk 
path must pass through (2j, 2fc) and in the next step intersect either ei or 62- It 
will go to cither {2j + 1, 2fc + 1) or (2j + 1. 2fc — 1) which means that wc obtain a new 
horizontal edge that is not covered by a dimer and we can repeat the argument. A 
random walk path connot inresect both ei and 62 because then they would have to 
meet at {2j, 2k) which would contradict the non-intersection condition. A similar 
argument applies at all locations in the graph. Hence, a dimer configuration gives 
rise to non-intersecting paths, and conversely if we have non-intersecting paths we 
can find the dimer covering corresponding to them. 

Consider a dimer covering of Gm,n- The total number of vertical dimers is 2ML, 
where L is the number of non-intersecting paths in the path description. The total 
number of dimers is M{2N -1-1) and hence the numbcer of horizontal dimers is 
M{2N -|- 1 — 2L). We will now define a probability measure on the set of dimer 
configurations by letting vertical dimers have weight w and horizontal dimers weight 
z. Let gM.N{infi,n) denote the number of configurations with m horizontal and to 
vertical dimers; each of these have the same probability. The partition function is 
given by 

(5.1) Z = Zm,n{z,w) = ^ 5(M,iv(m,n)^;"w". 

m , n > 

If Gi is the number of non- intersecting path configurations with exactly L paths, 
then we must also have 

N N 2MI, 

(5.2) Z=Y, GL^^(2^+l-2i)^„2ML ^ ^M(2iV+l) ^ 

L=0 L=0 ^ 

Recall that the possible initial (=final) positions for the non-intersecting paths are 
{2k \ < k < N}. Let Gl{x) denote the number of configurations with L non- 
intersecting paths whose initial conditions are 2x\ < ■ ■ ■ < 2xl, where xi, . . . ,xl S 
[N] = {0,...,N} are given. Then, 

(5.3) Gl= Y1 ^^(^)- 

0<xi<-<xz,<N 

Hence, the probability that the initial positions x given that we have exactly L 
non- intersecting paths is Gl{x)/Gl, and we define a probability on [N]^ by 

(5.4) ul{x) = j^Gl{x), 

where the right hand side is extended to [iV]^ by requiring it to be a symmet- 
ric function. The I -particle correlation function, given that the total number of 
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particles is i, is defined by 



(5.5) Ri,l{xi, ...,xe)= . ^' ^, ^ ul{xi, . . . ,xl)- 



xi^i,...,XLe[N\ 

The probability of having exactly L particles is, by ( |5.2| ) 

%(u;/z)2*"z*^(2Ar+l) 

and hence the (-particle correlation function, with no restriction on the total number 
of particles, is 

^M(2Af+l) N 

(5.6) Riixi,...,Xi) = ^i?,,L(a:i,...,x,)(zz;/z)2A^^GL. 



Z 

L=£ 



Set 



(5-7) (/-(s, i) = Y 

<s,t< N, where c™ ^ 1/2 if m = or to = 2Ar and = 1 if 1 < to < 2iV - 1. 
Proposition 5.1. Set 

N 



i2w/zy^Wj 

j=0 \ I J 3 



(5.8) K{x,y)=J2cl^{x,my,j)- 
where wj = (cos^)^*'^. Then, 

(5.9) . . . ,a;£) = det{K{x^, Xj))ij^T^. 
Also, 

(5.10) Zm,^(z,H = 11(1 + (^!^cos||)^*0. 

k=0 

Proof. We consider a simple random walk on {— 1, 0, . . . , 2A^ + 1}, where — 1 
and 2N + 1 are absorbing barriers. The transition matrix is given by „ = if 
< n < 2A^, P* „ = 1 if n = -1 or n = 2iV + 1, P*.„_i = 1/2 if 1 < n < 2iV + 1 
Pn,n+i = 1/2 if < n < 27V - 1 and „ = otherwise. Since we are only 
interested in random walk paths that are not absorbed we can concentrate on the 
submatrix P — (P^ „)o<m,n<2Af ■ We are interested in the probability of going from 
2x to 2y in 2M steps, which is given by (P^*^)2x,2y A computation, see shows 
that 

N 

J2iP^^'h.,23Hj,k)=^kCb{x,k) 
3=0 

and because of the orthogonality 

N 

^<f>{s,j)(f>{j,t) = Sst, 
3=0 

we obtain 

N 

(5.11) (^'*')2.,2, = ^w;j0(x,jXj,y), 

3=0 
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with Wj as above. 

Now, it follows from the Karlin/McGregor theorem that 

(5.12) GUx) = 22^^Met((p2M)^^^^^jL^.^^^ 

where the factor 2^*^^ comes from the fact that Gjjx) counts the number of con- 
figurations. Note that the right hand side of ( 5.12 ) is a symmetric function of 
xi, . . . ,xl- Let / be a given function on N and set 

L 

(5.13) $d/]= <i^iiiP''%^.,2.,)f,,=^l[f{xj), 

xe[N]L j=i 

and wr ite [A^ lf (A ) for t he co efficient of in the power series F{X). It follows from 
(U), dU), ( Ml ) and ( p^ that 



(5.14) 

Also, we set 
(5.15) 



m 



2ML 



N 



LIGl 



[A>l[1 + A/]. 



L=Q 



and note that by (fj), ( p^ and ( p3| ) 

Z 



(5.16) 



$[1] 



2M(27V+1) ■ 



Now, by (5.6), (5.14), (5.15) and (5.16), 



(5.17) 



Rdxu...,x,)f{f{x,)^l\[X']^^^±^. 



By (5.11) and a classical identity, see [ |67[ , 

det((p2M)2^^^2^jL^.^^ = det(^u;fe0(a;„/c)0(fc,a;,)),^^.^i 

fc=0 

1 ^ 
= Z! ^ [det(0(a:,,%),^^.^i]2 -Q^^^^ 

Thus, if we use the same identity again in the other direction we obtain 

N 

(5.18) $l[/]= det{YwkMx,k,)<f>{x,k,)f{x))f^^^,. 



Set 



x=0 



N 



t) ^Y '^«<^(^' t)f{x). 

Then, by ( |5.15D , ( ^.18D and a Frcdholm expansion 

(5.19) $[/] = det(5,,, + (2^«/z)2^'^/C(^, j))5=o- 
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Set ^ - ((1 + {2w/zf^"w,)5,,,)^^^^ and B = ((2u;/z)2*^/C(z, j))f^.^o. Then, by 

(iH), 



(5.20) 



$[1 + A/] =det(^ + AB). 



Combining ( 5.16 ) and ( 5.2C ) we obtain ( 5.10 ). Furthermore, ( 5.17 ) can be written 
(5.21) 



E 



Ri{xi, 



£\[\^] det(/ + XA-^B) = £l[X^] det(/ + XKf) 



with K given by (5.S). The last step is proved as in pTt, see also [p9l. Expanding the 



last expression in (5.21) in a Fredholm expansion gives (5.9) since / was arbitrary, 
and the proposition is proved □ 

We will now discuss the asymptotics of the partition func tion a nd the correlation 
functions. It is well known how to obtain a formula like ( 5.10| ) for the partition 
function using Kasteleyn's method, |Q. Also the two-point correlation function 
has been computed in |6^] using the methods of |l^ . See also 1 20 . The computation 
above is different and emphasizes the similarity between certain aspects of the dimer 
model and random matrix theory. The limiting correlation functions we obtain are 
given by the discrete or ordinary sine kernel. 

The number of vertices is 2M(2N +1) and the free energy per vertex is defined 

by 

1 



(5.22) 



fM,N{z,w) 



2M(2N + i; 



\ogZM,N{z,w). 



Proposition 5.2. The limiting free energy is 



fiz,w) 



lim lim fM.N{z,w) 

M^oo A'-.oo 



(5.23) 



- ; 2 



\\ogz 



ifw/z < 1/2 



i log z + 2^ log(2H cos s)ds if w/z > 1/2, 



Proof By ( |5.10| ) and ( |5.22| ) we obtain 



TV 

lim fMN{z,w) — -\oe-z^ hm > logil + ( — cos ) ) 

k=0 



- log z H 

2 ^ 2ttM 



""/^ 2w 

l0g(l + (— C0Ss)2A^)ds 



If ui/z < 1/2, then the integrand goes to uniformly and we obtain the first part of 
(5.22). If w/z > 1/2, then the integrand goes to unless < s < arccos(2;/2u') — 
t:9o/2 and we obtain our result by using the inequalities 



cos sf" < 1 + ( V - 2( V "^"^ 



(5.24) 

forO<s<6'o. □ 

If we only had horizontal dimers {L = 0), then Z = and hence the free 

energy per vertex is ^ logz. We can thus interpret the phase transition in ( 5.23 ) as 



saying that for w/z < 1/2 the system is completely frozen, whereas for w/z > 1/2 
we have both horizontal and vertical dimers. This type of phase transition is called 
a K-type transition in |53]. We will study the limiting correlation functions in the 
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non- frozen phase wjz > 1/2 and in a scaHng hmit where we approach the critical 
point from above. In this scaling limit we will obtain the sine kernel determinantal 
point process of random matrix theory. By proposition 5.1, formula (5.!;) it suffices 
to investigate the asyniptotics of the kernel (5.8). 



Proposition 5.3. Assume that wjz > 1/2. T/ien 



(5.25) 



N N 
lim lim K( hi, — 

Af^ooW^oo 2 2 



sin 7r(i — s)6q 



TT{t-s) ' 

for any fixed t,s ^ Z. Here Oq = ^ arccos ^ as before. Let cn, N > 1, be a given 
sequence such that e — *■ and N ^JTn — > oo as N oo. Assume furthermore that 
M = M{N) oo and Mcn oo as N ^ co. If2w/z = 1 + cat, then 

^2 ^2y/2^' 2 ^2y/2^'^ 

for any fixed ^, 77 G M. 



(5.26) 



7r(^ - v) 



Note that 6*0 gives the local density of the intersection of the random walk paths 
with a vertical axis and that — s- as w/z 1/2+. 

Proof. Assume that N ~2n for simplicity. Then, by the definition of K{x,y), 
lim K(n + t,n + s) 

j—0 y z An ' 

Split the last sum into two depending on whether j is even or odd and compute 
the limits of the Riemann sums obtained. This gives 



[sin T:tu sin usu + cos -Ktu cos ■nsu] 







TTU\2M 



1 + (2Hcos^)2A^ 



du. 



With 6*0 as defined above we see that the limit of this expression as M 00 is 



cos 7r(i — s)udu 



and we have proved (5.25). Write 7Ar = 1 

N/2 



(5.27) 



-T 



sin( 



T^t{2j - 1) 
N 



sin 7r(i — s)6q 
' Tr{t~s) ' 

Cat = 2 w/z. We have that 
,7rs(2j-l), ' 



^s{2j-l) 



2M 



I sin(- 



IN COS 



N 



2'Ktj 27131 

COS — — cos ■ 



N 



N 



JN COS ^ 
1 + 7Ar cos ^ 



1 + 7Af cos 

2M 



2N 



where the error is negligible for large N. Note that (1 + eat) cos > 1 if (ap- 
proximately) j < Ny/2eN /tt. Hence the summation in (5.27) can be restricted to 
1 < j < N^/2ei^/Tr and in the limit we arc considering the right hand side of ( ^.27 ) 
becomes 

'■V2/n 2 



TT 



cos{—^{^ - ri)u)du 



sin7r(^ — rj) 
tt{£, - 77) 



□ 
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